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I review the problem of motion for small bodies in General Relativity, with an emphasis on devel- 
oping a self-consistent treatment of the gravitational self-force. An analysis of the various derivations 
extant in the literature leads me to formulate an asymptotic expansion in which the metric is ex- 
panded while a representative worldline is held fixed; I discuss the utility of this expansion for both 
exact point particles and asymptotically small bodies, contrasting it with a regular expansion in 
which both the metric and the worldline are expanded. Based on these preliminary analyses, 1 
present a general method of deriving self-consistent equations of motion for arbitrarily structured 
(sufficiently compact) small bodies. My method utilizes two expansions: an inner expansion that 
keeps the size of the body fixed, and an outer expansion that lets the body shrink while holding its 
worldline fixed. By imposing the Lorenz gauge, I express the global solution to the Einstein equa- 
tion in the outer expansion in terms of an integral over a worldtube of small radius surrounding the 
body. Appropriate boundary data on the tube are determined from a local-in-space expansion in a 
buffer region where both the inner and outer expansions are valid. This bufi'er-region expansion also 
results in an expression for the self-force in terms of irreducible pieces of the metric perturbation on 
the worldline. Based on the global solution, these pieces of the perturbation can be written in terms 
of a tail integral over the body's past history. This approach can be applied at any order to obtain 
a self-consistent approximation that is valid on long timescales, both near and far from the small 
body. I conclude by discussing possible extensions of my method and comparing it to alternative 
approaches. 



PACS numbers: 04.20.-q, 04.25.-g, 04.25. Nx, 04.30.Db 

I. INTRODUCTION 

The problem of motion is of tremendous historical im- 
portance in General Relativity, both theoretically and 
experimentally. In conceiving of the theory, Einstein was 
fundamentally concerned with explaining the motion of 
bodies solely in terms of the geometric relationships be- 
tween them. And much of the observational evidence for 
General Relativity — e.g., the deflection of light by mas- 
sive objects, the post-Newtonian effects in solar-system 
dynamics, and the slow decay of binary pulsar orbits due 
to the emission of gravitational waves — is tied to analyses 
of motion. 

Despite the historical importance of the problem, the- 
oretical treatments of it have largely been confined to 
two limiting regimes: first, the Newtonian limit of weak 
gravity and slow motion, in which Newton's laws of mo- 
tion and relativistic corrections to them can be derived 
for widely separated bodies; and second, the point parti- 
cle limit, in which the geodesic equation and corrections 
to it can be derived for bodies of asymptotically small 
mass and size. Study of the Newtonian limit was pio- 
neered by Einstein, Infeld, and Hoffmann [1, 2] and is 
now fully developed in post-Newtonian theory (see the 
reviews [3, 4] and references therein). Study of the point 
particle limit is less well developed, and it has typically 
focused on proving that at leading order, a small body 
behaves as a test particle, moving on a geodesic of some 
background spacetime (see, e.g., Refs. [5-10]). 

The advent of gravitational wave detectors such as 
LIGO [11] and LISA [12] has rapidly broadened the scope 
of this research. Because these detectors have the poten- 



tial to accurately measure the dynamics of bodies in re- 
gions of very strong gravity, there is now a pressing need 
to go beyond either the post-Newtonian or test-particle 
approximations. For example, in binary systems, which 
are potentially important sources of gravitational waves, 
the two bodies will emit gravitational radiation, thereby 
lose energy, and slowly spiral into one another (see Ref. 
[13] for an overview of such systems). Once the two bod- 
ies arc very near one another, the post-Newtonian ap- 
proximation breaks down. 

If the two bodies are of comparable mass, then their 
motion during these final stages must be determined 
via a numerical integration of Einstein's equation (see 
Refs. [14, 15] for reviews of numerical approaches to this 
problem). However, if one of the bodies is much less 
massive than the other, then the entire inspiral can be 
treated analytically, rather than numerically, by utiliz- 
ing the point particle limit. (See Refs. [16, 17] for an 
overview of these systems, called extreme mass ratio in- 
spirals (EMRIs).) In this case, an expansion in the point 
particle limit roughly corresponds to an expansion of the 
metric in powers of the mass ratio m/M ^ e. At lead- 
ing order in this expansion, the small body moves on a 
geodesic of the spacetime of the large body, while simul- 
taneously emitting gravitational waves. Obviously this 
approximation breaks down after a brief time, since it 
implies that the body will travel forever on a geodesic 
even as it emits waves that carry off energy and angular 
momentum. Thus, one must proceed beyond the leading- 
order, geodesic approximation. At subleading order, the 
metric perturbations produced by the small body force it 
onto an accelerated worldline that slowly spirals into the 
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large body. The acceleration of this worldline, caused by 
the body's interaction with its own gravitational field, is 
called the gravitational self-force. Along with all other 
corrections to the test-mass approximation, this force will 
be the principal subject of the present paper. 

A formal expression for the gravitational self-force was 
first derived in 1996 by Mino, Sasaki, and Tanaka [18] 
and Quinn and Wald [19]; the resulting equation of mo- 
tion is now known as the MiSaTaQuWa equation. Since 
then, numerous other derivations have been offered (sec 
Refs. [20-27] and the reviews [28, 29] for examples). How- 
ever, none amongst this plethora of derivations has over- 
come a fundamental difficulty in defining a point particle 
limit: how can one accurately, self-consistently, and sys- 
tematically incorporate corrections into a worldline? 

At its most fundamental level, the self-force problem 
consists of finding a pair (7, /i^i/) representing the world- 
line 7 and metric perturbation of an asymptotically 
small body. This problem is far from trivial. Unlike in 
other field theories such as electrodynamics, equations of 
motion arc not independent of the field equation in Gen- 
eral Relativity — in fact, they are integrability conditions 
for the Einstein field equation, following from the restric- 
tion imposed by the Bianchi identity [1, 28]. This means 
that at each order in perturbation theory, the equation of 
motion, and hence the worldline, is fixed by the Bianchi 
identity; using any other worldline means that a given 
nth-order perturbation /i'^^ is not a solution to the nth- 
order Einstein equation. But at each order, the worldline 
determined by the Bianchi identity differs from that at 
every other order. It seems clear that the higher-order 
equations of motion are corrections to the lower-order 
ones, but there is no obvious way to self-consistently and 
systematically incorporate these corrections. 

For example, suppose that in an EMRI system, one 
expands the metric and the Einstein equation in pow- 
ers of e, and that at order e the stress-energy tensor 
of the small body can be approximated by that of a 
point mass. Then the linearized Einstein equation reads 
(5G^i/[/i] ~ 87rT^i^[7], where 5G^jv[/i] is linear in the per- 
turbation /la/3 and r^i/ [7] is the stress-energy tensor of a 
point mass moving on a worldline 7 in the background 
spacetime. The linearized Bianchi identity implies that 
the point particle source must be conserved, which in 
turn implies that the particle must move on a geodesic 
of the background spacetime. Thus, at first order in the 
expansion, the body travels forever on a geodesic. In or- 
der to correctly derive the self-force, we must proceed to 
second order. But even if that is accomplished, it seems 
that the self-force cannot trivially be identified as the 
equation of motion of the worldline, since the worldline 
is fixed by the first-order Bianchi identity. 

As in any problem involving a small parameter, two 
options present themselves: first, one can assume a regu- 
lar Taylor series expansion of every function in the prob- 
lem, which leads to a succession of equations that can 
be solved exactly, order- by-order; or second, one can be 
satisfied with the construction of an approximate solu- 



tion to the exact equation, however the solution may be 
arrived at. If the first approach is adopted, then the 
linearized Bianchi identity fixes the worldline to be a 
geodesic, and the self-force can only be interpreted in a 
perturbative sense, as the equation of motion for a "devi- 
ation vector" describing an infinitesimal correction to the 
geodesic. But this interpretation is meaningful only for a 
brief time: since the particle will eventually plunge into 
the large body, the "small" correction to the geodesic will 
eventually grow large. At that point, the entire expansion 
in powers of e will have broken down. In other words, the 
straightforward expansion of the Einstein equation is a 
valid approximation to the actual Einstein equation only 
on short timescales. 

However, in studies of the problem of motion in Gen- 
eral Relativity, this first approach is atypical; instead, 
the second approach is typically adopted. In the self- 
force problem, this has been realized in the procedure 
of "gauge relaxation," [19] in which the linearized Ein- 
stein equation is written in the Lorenz gauge, leading to 
a wave equation that can be solved for a point particle 
source moving on an arbitrary worldline — but then, in 
order to circumvent the conclusion that that worldline 
must be a geodesic, the solution to the linear equation 
is allowed to slightly violate the Lorenz gauge condition 
[18, 19]. This procedure yields an approximate solution 
to the exact Einstein equation, as desired, and it leads 
to a single worldline obeying a self-consistent equation of 
motion. It is also similar to successful methods of post- 
Newtonian theory, in which, prior to any expansion, the 
exact Einstein equation is written in a "relaxed" form 
that can be solved for an arbitrary source. However, 
the gauge-relaxation used in the self-force problem lacks 
the systematic nature of the post-Newtonian method, in 
the sense that the relaxed linear equation has not been 
shown to follow from a systematic expansion of the exact 
Einstein equation, and the solution to the relaxed linear 
problem has not been related to a solution to the exact 
problem. 

One of the goals of this paper is to provide such a sys- 
tematic justification of the gauge-relaxation procedure, 
by breaking the exact Einstein equation into a sequence 
of equations that can be solved for a fixed, e-dependent 
worldline. Of course, trying to do so introduces another 
problem: a point particle is a sensible source only for the 
first-order, linearized Einstein equation. In the full, non- 
linear theory, point particles are not mathematically well- 
defined sources.^ Going beyond first order thus means 
that one must abandon the point particle approximation. 

Hence, we can see that an accurate, self-consistent so- 



At least this is true within classical distribution theory [30] , since 
the Einstein tensor of a point particle would contain products 
of delta distributions and hence be too singular to be treated 
as a distribution. However, more general methods based on 
Colombeau algebras, which allow for multiplication of distribu- 
tions, have been devised to overcome this problem [31]. 
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lution to the self-force problem should accomplish all of 
the following: 

1. go beyond first order pertm'bation theory, 

2. define a self-consistent, "corrected" worldlinc, 

3. treat the body as extended, but asymptotically 
small. 

Note that these criteria are interdependent, since 
defining a meaningful worldline — one that incorporates 
corrections — requires going beyond first order, and going 
beyond first order requires accounting for the extension 
of the body. Furthermore, note that these theoretical 
goals are closely related to experimental ones: in order 
to extract the parameters of an EMRI from a gravita- 
tional wave signal, one requires a template that relates 
the signal to the motion of the body over a large portion 
of the inspiral. Such a template must be based on an ap- 
proximation scheme that is uniform on a domain of size 
~ in other words, the errors in the approximation 
must remain small over a long span of time. This can 
be accomplished only in a scheme that self-consistently 
incorporates the corrected motion. 

A. Organization of this paper 

The paper contains two main parts. The first part, 
comprising Sees. II and III, consists of a more thorough 
explication of the problem. This explication serves two 
purposes: first, to review the foundations of the problem 
and the various derivations in the literature. Much of this 
review overlaps with previous discussions by Mino [32- 
34], Hinderer and Flanagan [35], and Gralla and Wald 
[20]. However, because of the wealth of derivations per- 
formed to date, a more comprehensive review is timely. 
In addition, my presentation differs significantly from 
those earlier discussions and serves to motivate and pro- 
vide the necessary context for my own approach. The 
second purpose of the explication is to introduce the no- 
tion of a self-consistent expansion in which the metric 
perturbation is first written as a functional of a world- 
line and then expanded while holding the worldline fixed; 
in this expansion, the solution to the Einstein equation 
is consistent with what would result from an evolution in 
time that began with (1) some arbitrary initial data and 
(2) a system of evolution equations that involve only local 
values of the position, momentum, and metric perturba- 
tion of the small body at each value of time. Section II 
presents this expansion in the context of a point parti- 
cle; Section III then generalizes it to asymptotically small 
bodies. After laying that groundwork, the first part of 
the paper concludes in Sec. Ill C with an outline of my 
approach. Readers who are uninterested in the prelimi- 
nary material can skip directly to that section. 

The second part of the paper, comprising Sees. IV 
and V, presents my calculation of the first-order gravita- 
tional self-force and the metric perturbation created by 



an asymptotically small (sufficiently compact) massive 
object. Section IV presents the derivation of the first- 
order gravitational self-force; the result of this calculation 
is that the body moves on a geodesic of the spacetime 
g + , where is a homogenous perturbation that is 
regular on the body's worldline. In Sec. V, I calculate the 
metric perturbation induced by the body, which deter- 
mines in terms of tail integrals and recovers the usual 
MiSaTaQuWa equation. Section VI concludes the paper 
with a comparison to other methods and a discussion of 
higher-order and globally accurate approximations. 

Throughout this paper, Greek indices a,/?, 7,... run 
from to 3 and refer to a coordinate basis, uppercase 
Latin indices from the middle of the alphabet (/, J, K, ...) 
run from to 3 and refer to an orthonormal tetrad basis, 
lowercase Latin indices i,j,k,... run from 1 to 3 and refer 
to the spatial part of both the coordinate basis and the 
tetrad basis, and uppercase Latin indices from the be- 
ginning of the alphabet {A, B, C, ...) run from 1 to 2 and 
refer to angular coordinates. Sans-serif symbols g, R, T, ... 
denote exact quantities to be expanded, and denotes 
the covariant derivative compatible with an exact metric 
g; a semi-colon and V are used interchangeably to indi- 
cate a covariant derivative compatible with a background 
metric g. I work in geometrical units in which G = c = 1. 
I will frequently omit indices for simplicity. 



II. PRELIMINARY ANALYSIS I: THE MOTION 
OF A POINT PARTICLE 

Assume for the moment that the exact Einstein equa- 
tion, G^i/ = SttT^i, , can be made sense of with the point 
particle source 

T^^'lg,^]^ f mu''u''6{x,z{t))dt, (1) 
J J 

where 7 is the worldline of the particle, ^"(t) are the co- 
ordinates on 7, is its four-velocity, t is proper time 
with respect to g on 7, and S{x, x') = 5'^{x'^ — x'^)/ yj\g\ 
is a covariant delta function in the spacetime of the ex- 
act solution g^^, with |g| denoting the absolute value of 
the determinant of g^j/. The motion of the particle is 
constrained by the Bianchi identity ^^G'^'^ = 0, which 
implies the conservation equation ^^T^^'^ = 0. The con- 
servation equation in turn implies = {^^u^^) u'^ = 0; 
that is, 7 must be a geodesic in the full spacetime de- 
scribed by g. 



A. Non-systematic expansion 

Now suppose that m is small compared to all other 
length scales of the system, denoted collectively by TZ, 
and that we wish to construct an approximation to g^^^ 
and 7 for e = m/ inf 7?. <C 1. In this limit, the metric can 
be expanded as g = g -f- ft-, where h ^ e. And the exact 
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equation of motion a'^ = can then be expanded as [29] 

= -h^g""" +u^u''){2Kp,,~hp„,,)uPu^ + 0{e^) , (2) 

where a'^ = u^'-iyu'^ is the acceleration in the background 
spacetime g. Since the metric perturbation of a point 
particle will diverge as 1/r near 7, this equation of mo- 
tion is ill-behaved. Hence, some form of regularization is 
required. 

The very first derivation of the self-force, given by 
Mino, Sasaki, and Tanaka [18], followed earlier deriva- 
tions of the electromagnetic self-force [36, 37] by using 
conservation of energy-momentum, calculating the flux 
through a surface around the body and setting it equal 
to the change of energy-momentum within the tube. Un- 
fortunately, the regularization in this method essentially 
consists of discarding various divergent integrals and as- 
suming results for others [29]. 

Other derivations of the self-force [19, 22, 29] began 
by assuming that Eq. (2) is essentially vahd, but that 
only a certain regular part of h actually contributes 
to the acceleration. This regular part was assumed to 
be either the angle-averaged field [19] or the Detweiler- 
Whiting [22] regular field; the two assumptions yield the 
same force. All of these derivations are axiomatic — in 
the sense that they simply assume a form for the force — 
with their axioms supported by various plausibility argu- 
ments. Gal'tsov et al. [26] later showed that Eq. (2) can 
be regularized via a straightforward expansion along the 
worldline, without making any assumption about which 
part of the field contributes to the force. 

These derivations are based on solving the linearized 
Einstein equation, substituting it into Eq. (2), and then 
regularizing the result in one way or another. Expanding 
the exact Einstein equation to linear order in e, we find 
at zeroth order that G^'' = 0, which tells us that the 
background metric g is that of a vacuum. At linear order 
we find 

6G^"'[h]^87:T^"'[-f], (3) 

where SG'^'^ is the linearized Einstein tensor, and r^'^[7] 
is the stress-energy tensor of a point-mass moving on a 
worldline 7 in the background spacetime defined by g. A 
formal solution to this equation can easily be obtained 
by imposing the Lorenz gauge condition. In the Lorenz 
gauge, the first-order Einstein equation is split into a 
wave equation, which I will write as 

E^,[h] ^ -IGttTM (4) 

and the gauge condition '^^Ti^^y = 0, which I will write 
as 

L^[h] = 0, (5) 
where iJ^j^ and are linear operators defined by 

EpM = «5^V^V^ + 2i?//) V, (6) 



An overbar indicates trace-reversal with respect to g; e.g., 
h = h — io oP'^h 

The retarded solution to the wave equation is 

= 4my"G^,^V'u'''u"'rfi', (8) 

7 

where G^^f^iyi is the retarded Green's function for the 
operator J?^^ and t' is the proper time with respect to g 
on the worldline. (My conventions for Green's functions, 
along with useful identities, are given in Appendix A.) 
Near the worldline, this solution can be decomposed into 
a local term, which diverges as 1/r, and a so-called tail 

tail 

term h , defined as 

tail /■*■■<=' - ' ' , 

hp^ = 4m / G^^^v'W^ u"" dt', (9) 

where the upper limit of the integral is cut off just prior 
to the retarded time t^ct in order to avoid the divergence 
of the Green's function there. The regularized equation 
of motion is obtained by replacing the exact field h with 

tail 

the regular field h in Eq. (2). 

But in order for the solution to the wave equation to 
also be a solution to the first-order Einstein equation, it 
must satisfy the gauge condition, which now reads 

= L^h] 
= 4y V'^G^,^v'T'^'"'[7]rf^' 

= ^J G^^,%,T'^''''[j]dr 

= 4m f G^^^^a^'dt', (10) 

where I have used the identity V^G^i/^'jy' = —G^(^^'-^i) 
(Eq. (A6)) and integrated by parts in going from the sec- 
ond line to the third. Thus, we find that imposing the 
gauge condition is equivalent to imposing the conserva- 
tion equation T^^"^.^ = 0, which is equivalent to imposing 
the first-order Bianchi identity SG^'^-i, = 0. (The same 
equivalences could also be found by taking the divergence 
of the wave equation.) The consequence of any of these 
conditions is that 7 must be a geodesic in the background 
spacetime. 

This requirement obviously contradicts the equation of 
motion (2), regularized or not. In the earliest derivations 
of the self-force [18, 19], the contradiction was overcome 
by allowing the Lorenz gauge to be slightly violated by 
the first-order perturbation, effectively sidestepping the 
requirement that 7 must be a geodesic. This would mean 
that the metric perturbation ft, is a solution to the wave 
equation (4) but not a solution to the linearized Einstein 
equation. Besides the fact that this gauge-relaxation pre- 
vents h from exactly solving the linearized field equation. 
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it also calls into question the expression for the self-force. 
The force is presumably an integrability condition for the 
second-order Einstein equation, but the point particle 
source is meaningful only at linear order. By taking the 
accelerated motion into account, we implicitly assume 
the solubility of the second-order problem, even though 
a point particle source ceases to makes sense at second 
order. 

As Mino, Sasaki, and Tanaka stated in their original 
derivation [18], neither of these problems are too severe. 
The failure to exactly solve the first-order Einstein equa- 
tion is not devastating, since if the acceleration is of or- 
der e, then the errors in the gauge condition are of order 
presumably these errors would be cancelled at second 
order. And we might comfortably presume that if the 
extension of the body were somehow taken into account 
at second order, the equation of motion would remain 
consistent with the MiSaTaQuWa result, which is, after 
all, fairly well motivated. However, can we justify these 
reassurances more systematically? In the remainder of 
this section, I will examine methods of self-consistently 
incorporating the acceleration into the point particle so- 
lution. In the next section, I will begin to discuss how 
these methods can be applied to an asymptotically small 
extended body. 

But first, allow me to introduce some nomenclature. In 
general, for any function f{x,e), there are two types of 
expansions to consider: a regular expansion, of the form 

N 

/(x,e) = 5]£"f"'(-^)+0(£'^+'), (11) 

n=0 

where the cocfEcients /<"'(a;) are independent of and 
a singular expansion, of the form 

N 

/(x, e)-^ (a:, £) + 0(e^+i), (12) 

n=0 

where the coefficients f'-"\x,e) depend on e but are of 
order 1, in the sense that there exist positive constants k 
and Eq such that [/'"'(a:, £)| < A: for < e < eo, but 
lim /'"'(a;, e) ^ (unless /<">(a;,e) is itself identically 

e— >0 

zero). Put simply, the goal of a singular expansion is 
to expand only part of a function's £-dependence, while 
holding fixed some specific e-dependence that captures 
one or more of the function's essential features. Further 
details of such expansions, and their role in singular per- 
turbation theory, can be found in numerous textbooks 
[38-41]; see the text by Eckhaus [42] for a rigorous treat- 
ment of some aspects. See Rcfs. [43, 44] for discussions 
of singular perturbation theory in the context of General 
Relativity. 



^ More generally, a regular expansion can be of the form f{x,£) = 
X]^JLo -'^i(^)/'"' (^) + 0(\n+i) for any sot of functions An sat- 
isfying limg^o An+i(£)/An(e) = 0. An analogous generalization 
holds for singular expansions. 



B. Regular expansion 

Let us begin by considering a regular expansion in pow- 
ers of e. Although such an expansion might be at the 
backs of most researchers' minds when they derive an 
expression for the self-force, only one extant derivation 
of the MiSaTaQuWa equation [20] has explicitly sought 
to remain within the framework of such an expansion. 
We begin by expanding the metric as 

gt,y{x,e) = g{x) + ehllix) + e^ll'J^ix) + 0{e^) , (13) 
and the Einstein tensor as 

+ e^S^G'"'[h^''] + 0{e'^) , (14) 

where G is the Einstein tensor of the background metric 
g, SG^'^[h] is linear in h and its derivatives, and S'^G^'^[h] 
is quadratic in them. Similarly, by expanding the \/[g[ 
that appears in it, and converting from the proper time 
in g to the proper time in g, the stress-energy tensor can 
be expanded as 

T'^''[g,7] = eT^''[-/]+e^6T^^''[}^"\j] + 0{e') , (15) 

where the factor of e is pulled out of T for clarity. How- 
ever, given that 7 satisfies the geodesic equation a*^ = 
in the full spacetimc, it should be obvious that it will 
generically depend on e, so the above expansion is not 
yet regular. To make it regular, we must expand the 
worldline as 7 = 7'"' -I- £7'^' -I- O(e^). With the coordi- 
nates of the worldline defined by z^lt, e), this expansion 
takes the form 

z'^it,s) = z^,^it) + sz^,^{t) + 0{s^), (16) 

where for the remainder of this subsection, t will indi- 
cate proper time on the leading-order worldline 7"". To 
make this expansion most meaningful, we can insist that 
at some time t = to the exact curve z" is tangential to 
the leading-order curve z"gy the corrections z"j, n > 1, 
then determine the deviation of the exact curve from the 
geodesic as time progresses away from t = to- Since the 
different terms in the expansion cannot map to different 
points in a curved spacetime, the "corrections" are in fact 
vectors defined on the leading-order worldline; they "con- 
nect" Z("j (i) to z" (t) , in the same sense as a geodesic devi- 
ation vector connects two neighbouring geodesies. Hence, 
in this expansion, one does not arrive at an equation for 
the acceleration of a worldline. Instead, one arrives at 
an equation for the acceleration of a deviation vector. 
This acceleration will naturally include a term identical 
to that of the geodesic deviation equation [20] , due to the 
drift of the the true worldline 7 away from the reference 
worldline 7'"'. 

By using this expansion of the worldline, we can con- 
struct a regular expansion of the stress-energy tensor, 

+ e'~6T,4r\l''']+0{e'), (17) 
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where ST'^" is linear in and ST^" is linear in where z 

Substituting this expansion into the Einstein equa- 
tion, we arrive at a sequence of field equations, written 
schematically as 



(1) 



G"" = 0, 



(18) 
(19) 



These equations can be solved order-by-order for the 
background metric g, the perturbations h^J^, and the 
curves z^„y 

At linear order, repeating the analysis above, we can 
split the Einstein equation into a wave equation and a 
gauge condition, and the first-order perturbation li ' can 
be written as 



, , ii' j-i' 



(21) 



where z^a^ = ^^"^ . Imposition of either the gauge con- 
dition, the Bianchi identity, or the conservation of the 
source implies that 7'°' must be a geodesic in the back- 
ground spacetimc. 

But 7*°' does not describe the true worldline of the par- 
ticle: the effect of radiation at first order is incorporated 
into the correction zl{^ , which is again determined by the 
Bianchi identity (or, equivalently, the conservation of the 
source 87r(5T^''-h87r(5r^''-(52G^''). Explicitly, after split- 
ting the second-order field equation into a wave equation 
and a gauge condition, the solution to the wave equation 
is given by 



flu 



G 



7/7/(7/ 7/ -o'''"'\h , ,cH' 



r , , (0 / 7'^' 



. ' fi ' u' • (5' • 7 

) + ^(o)^(o)57'<5' ^(0)^(1) 



^ , , J.I' „P Ji/ 

flVfJ,' u' -.p' ^(0) ^(0) '^(1) 



:f / Gp^^^.^S^G^^" dV, 



• u 

%) = 



The first line in this solution arises 
from ST, while the second and third arise from ST. 

(2) 

Imposing the gauge condition L|^^[h ] = 0, making use 
of Eq. (A6), integrating by parts, and then making use of 
the Ricci identity and the second-order Bianchi identity 
(given by V^^^G^" = -^E^^^G^t - S^^p^SG^''>, where 
6r is the linear correction to the background Christoffel 
symbol) we arrive at 



^(0)^(0) -^(i)- "T^^^ equation describes the 
spatial deviation of the true worldline away from the 
reference geodesic because it is an equation for a 

deviation vector, it includes a term proportional to the 
Riemann tensor. 

While the basic idea of this approach is valid and rig- 
orous, it is unsatisfactory because of its limited realm of 
validity. For example, in a typical EMRI orbit, the radial 
coordinate z'^^^ on the particle's leading-order worldline 
(20) will be of order e° for all time, while the deviation vec- 
tor z^jj will grow as a'"^' • {t — to). This means that the 
expansion of 7 is valid only on timescales t ^ TZ: after 
a radiation-reaction time t ~ TZ/s, the "correction" £Z(i) 
will be of the same order as the leading-order term Zj,,) . In 
other words, the expansion is not uniform in time. And 
once we commit ourselves to a nonuniform expansion, 
we must restrict the entire problem to a bounded time- 
interval [ti,tf]. Within this fixed interval, the expan- 
sion is valid in the sense that we can guarantee that our 
approximation will be accurate to any given numerical 
value by making e sufficiently small; on an unbounded, or 
a generically e-dependcnt interval, this statement would 
not hold true. The restriction to this bounded region 
has several important consequences. Most obviously, as 
previously stated, we are specifically interested in large 
changes that occur on the time-interval ~ TZ/e — such as 
the particle's slow inspiral in an EMRI. Thus, the entire 
expansion scheme fails on the timescale of interest. 

The restriction to a bounded time-interval also re- 
stricts the formalism in an important way: since the 
expansion of the Einstein equation is valid only on a 
bounded region, the solution to it cannot (necessarily) be 
expressed in terms of an unbounded past history. In other 
words, the boundedness of the domain effectively forces 
us to cast the problem in an initial value formulation from 
the beginning. This means that we cannot express the 
force purely in terms of the usual tail integral; as soon as 
one writes down the solution as an integral over the entire 
past history, one assumes that one's expansion is globally 
valid, rather than just locally valid. ^ We can easily see 
this from the following argument: The correction terms 
Z(„) grow large not only for times far in the future of tg, 
but also for times far in the past of to- Hence, at any time 
t, the difference between the tail as calculated on 7*°' and 
the tail as calculated on 7'°' +eY^^ will differ by a signif- 



dt' 



(22) 



I'pcr -^(0)^(0)^(1) 



9(5^ + ^('o)^(o))(2^i/p: 



t;j')-^(0)^(0)i (23) 



icant amount, given by 

z^^^{t')da'G{x, Zi^o){t'))dt'; since Z(i) grows with 
\t — to\, the difference between the two tails appears to 
be potentially infinite. It is quite likely that the decay of 
the retarded Green's function would ameliorate this di- 



^ I have assumed for simplicity that 2:^^^ is perpendicular to i^^^y 
In the general case, the left-hand side of the equation reads 
{9'^" ~^ ^{o)^(o)')9'^p'^{i)' siich that the result still describes only 
the spatial part of the acceleration. 

* This point seems to have been missed in Ref. [20]. 
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vergence in any case of interest. But there is no obvious 
reason for extending the domain of the solution beyond 
the domain of vahdity of the expansion. 

Hence, at each order, the integral over the source must 
be cut off at the initial time t = ti, and the remainder of 
the tail must be replaced by Cauchy data on that initial 
timeslice. A consequence of this is that the self-force is 
not naturally expressed in terms of a tail integral over 
an infinite past history; instead, it is more naturally ex- 
pressed in terms of a purely local regular field, defined 
as the retarded field minus a certain local, singular part, 
in the manner of Detweiler and Whiting. Besides mak- 
ing the solution valid, this also has the advantage of ex- 
pressing the force in terms of local quantities, with no 
reference to the past history of the particle; this is use- 
ful for a numerical integration in the time domain — and 
for developing a two-timescalc method, as I will discuss 
presently. 

C. Singular expansion 

Given the limitations of a regular expansion, let us now 
consider a singular expansion. In effect, this expansion 
will formalize the non-systematic procedure presented in 
Sec. II A and provide a systematic justification of the 
relaxation of the Lorenz gauge in the first-order problem. 

Recall that our basic goal is to find a pair (7, /i) sat- 
isfying Einstein's equation. In a regular expansion, both 
the worldline and the metric perturbation are expanded 
in the limit of small e. In the singular expansion we shall 
now consider, the worldline is held fixed. To find each 
term in this expansion, I seek an expansion of both the 
exact Einstein equation and the exact equation of mo- 
tion that can be solved with this fixed worldline. Hence, 
I decompose the metric as 

g,,^(x,e) + /i^^(x,£;7), (24) 

where the semicolon is used to separate ordinary 
coordinate-dependence from functional dependence; 
when leaving the dependence on coordinates implicit, I 
will write, e.g., h^v[')]. I assume that the perturbation 
can be expanded while holding this functional depen- 
dence fixed: 

N 

K.{x-n) - ^e"C(^;7) + 0(£^+') , (25) 

n=l 

where each term li^^^ is a functional of the true world- 
line 7 but is nevertheless of order unity. Note that the 
approximation scheme fails — becoming both inaccurate 
and internally inconsistent — if any of these coefficients 
are found to grow larger than order unity. Substituting 
Eq. (24) into the Einstein equation, we arrive at 

5G^"' [h] =-G'"' + SneT'"' [7] + SneST'''' [h, 7] 

-6^G^"'[h]+ ... (26) 



Now, one might seek to expand and solve this equa- 
tion at fixed 7. Unfortunately, if one substitutes the 
expansion (25) into Eq. (26), then one arrives at the 
linear equation (5G^''[/i'*] = 87rr^''[7]. And since the 

linearized Bianchi identity holds for any h , this equa- 
tion implies that '%T^'^[j] = — and hence that 7 must 
be a geodesic. Thus, this method immediately fails. To 
avoid such a problem, I recast the full equation in a more 
useful form by imposing the Lorenz gauge condition on 
the entire perturbation h (rather than on any individual 
term in its expansion):^ 

L^[h]=0. (27) 

This transforms the Einstein equation into the weakly 
nonlinear wave equation 

E^"'[h] = 2G^'' - 167r£r'^'^[7] - 16^^eST^"'[h,-/] 

+ 25^G''''[h]+ ... (28) 

Unlike Eq. (26), this equation can be solved for an ar- 
bitrary worldline. Its form is essentially identical to the 
relaxed Einstein equation, which forms the basis of most 
post-Minkowski expansions [45, 46]. Both equations are 
"relaxed" in the sense that they can be solved without 
specifiying the motion of the source; the motion of the 
source is determined by the gauge condition. Also, in 
both cases, nonlinearities are treated as source terms for 
a hyperbolic wave operator; this means that corrections 
to the null cones are incorporated into the perturbations, 
rather than into the characteristics of the wave equation. 

Substituting the expansion (25) into the wave-like 
equation (28) and solving at fixed 7 now yields the se- 
quence of equations 

Gf"" = 0, (29) 
E^'•'[Jl'^] = -IGTrT^^'i-f], (30) 
£:^''[/i'"] = -16TTdT'"'[h''\j] +2d^G^"'[h^''] (31) 

Along with this expansion of the Einstein equation, I 
seek an analogous expansion of the matter equation of 
motion = 0: 

(i , e) = ajT' (t) + ea)i' (i ; 7) + O (e') ■ (32) 



I assume that this condition can be imposed up to any desired 
order of accuracy. At second order, it can be imposed via a 
gauge transformation generated by a gauge vector f satisfying 
= L°'[h] + S"[h,S.] + 0{e^), where is quadratic in h 
and ^. At higher order, higher-order terms would appear on the 
right-hand side. Beginning with h in an arbitrary gauge, this 
weakly nonlinear wave equation can be solved iteratively for ^, 
using the same methods as those used to solve for h. Note that 
in the singular expansion presented here, gauge transformations 
have the coordinate form — > — e^°'{x,'y) + O(e^). 
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This is an expansion of a function of time along the fixed 
worldhne. Of course, it must also follow from the Bianchi 
identity, and it must therefore also follow from the Lorenz 
gauge condition. Hence, I seek an expansion of the gauge 
condition that is equivalent to the expansion of the ac- 
celeration. Substituting the expansions of the perturba- 
tion and the acceleration into the exact gauge condition 
Lfj^[h] = 0, and solving at fixed 7, yields the sequence of 
equations 



e: 



L<"'[r]=o, 



(33) 
(34) 



where = L[f]\^^^,,„ L^'^[f] is linear in a<^', L<^'[/] 

is linear in a'^' and quadratic in a*^\ and so on. Note the 
meaning of the notation here: Ii ' [7] is the coefficient of 
at fixed 7 — including, in particular, at fixed values 
of the acceleration on 7 — while for / ~ 1, is the 

coefficient of e" in an expansion that incorporates an 
expansion of a on 7. 

This expansion of the gauge condition serves to de- 
termine increasingly accurate equations of motion for 
the single fixed worldline 7. Alternatively, it might be 
thought of as an iterative improvement of the choice of 
worldline. In cither case, in any concrete calculation, the 
metric perturbations would be treated as functionals of 
the worldline described by the highest-order equation of 
motion available. 

The solution to the first-order wave equation, (30), is 
given by 



(35) 



where t' and are, respectively, the proper time and 
four-velocity on 7. Note that this is the "usual" solu- 
tion obtained by solving the linearized wave equation, as 
in Sec. II A. The acceleration of the true worldline is 
determined from Eq. (33): 



= L<»'[/^'1 



4m / G^'^ a™dt, 



(36) 



which implies that a*J'' = 0. It docs not imply that = 
0, however. 

Proceeding to second order, the solution to Eq. (31) is 
CW=2 / G^.^vm'^ «'''(«'' u^'-/^>';',,dt' 



^ / G^.^v'^'G^'^^'dT/'. 



(37) 



Imposing the gauge condition (34), making use of 
Eq. (A6) and the second-order Bianchi identity, and in- 
tegrating by parts determines the acceleration to order 



a=0 



(38) 



Note that the right-hand side of this equation is eval- 
uated on the worldline, and once evaluated, it contains 
a term proportional to — o./^, corresponding to the anti- 
damping phenomenon discovered by Havas [47] (as cor- 
rected by Havas and Goldberg [48]). However, my as- 
sumed expansion of the acceleration has forced the right- 
hand side to be evaluated for a ~ a'"' = 0, which serves 
to automatically yield an "order-reduced" equation with 
no higher-order derivatives (assuming, of course, that a 
time-derivative does not change the order of a term in 
the expansion). 

From this calculation, we see that this method yields 
an equation of motion that agrees with the expansion 
given in Eq. (2). In both cases, the equation of motion 
applies to the actual worldline 7, not to a correction to a 
reference geodesic. Combined with the first-order pertur- 
bation given in Eq. (35), the equation of motion defines 
a self-consistent solution to the Einstein equation, up to 
errors of order on a timescale TZ/e; combined with the 
sum of the first- and second-order perturbations, it de- 
fines a solution accurate up to errors of order on a 
timescale of order TZ. 

One should note two important facts about the results 
just derived. First, from these results, one can easily 
derive those of the regular expansion, given in Eqs. (23) 
and (21), by expanding the worldline and following the 
usual steps involved in deriving the geodesic deviation 
equation. Second, while the Lorenz gauge is especially 
useful for finding the metric perturbation in the singular 
expansion, it is not essential for finding the equation of 
motion, which could have equivalently been found from 
the conservation of the source 87re(r^'' + ST^"') - S^C". 

Beyond these specifics, one should also note the broad 
similarity between this singular expansion and a post- 
Minkowksian expansion (in particular, the fast-motion 
approximation [45]): the split of the Einstein equation 
into a wave equation and a gauge condition, the iterative 
solution to the wave equation in terms of an arbitrary 
worldline, and use of the gauge condition to fix the world- 
line. Given these commonalities and the many successes 
of the post-Minkowskian expansion, one might hope that 
the singular expansion suggested here will be equally suc- 
cessful in more general contexts. Note, however, that the 
character of the solutions given in Eqs. (35) and (37) is 
significantly different in a curved background than in a 
fiat one, since curvature creates caustics in null cones and 
allows gravitational perturbations to propagate within, 
not just on, those cones. These complications suggest 
that the integrals in Eqs. (35) and (37) might much more 
easily display secular growth in a curved spacetimc. Since 
the expansion is consistent only in the absence of such 
secular behaviour, it may be valid only in certain space- 
times and with certain initial conditions. 

Perhaps a more significant difference between the 
above expansion and a post-Minkowskian one is the 
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choice of gauge. The harmonic gauge used in post- 
Minkowksi expansions can be imposed as an exact coor- 
dinate condition ^Dx" = on the manifold of the exact 
solution; as long as the exact solution admits these coor- 
dinates, the gauge condition S^/i^'' = is automatically 
imposed on the entire metric perturbation, rather than 
on any particular order in its expansion. The Lorenz 
gauge used here, on the other hand, can be imposed only 
after decomposing the metric into a background plus per- 
turbation, and it is typically formulated only in terms of 
a first-order perturbation; there is no proof that it can 
be imposed on the total perturbation to all orders (al- 
though it can seemingly be imposed to any desired order 
by solving a weakly nonlinear wave equation). Despite 
these caveats, and independent of the analogy with post- 
Minkowskian theory, the expansion discussed here has 
the concrete advantage of offering a systematic justifica- 
tion of the self-consistent solution (35) and higher order 
corrections to it. 

Other arguments have been made in favor of using the 
self-consistent solution (35) rather than the regular solu- 
tion (21). The simplest argument is one based on "adi- 
abaticity": because the acceleration is very small, the 
true worldline deviates only very slowly from a geodesic, 
so the self-consistent solution can be "patched together" 
from a collection of regular solutions. This argument 
has been made frequently in the past, most recently by 
Gralla and Wald [20]. While it is intuitively reasonable, 
one must keep in mind its most basic assumption, which 
is that the (covariant derivative of the) tail integral as 
calculated over a geodesic 7*°' is nearly identical to the 
tail integral as calculated over the true worldline 7. This 
assumption is a very strong one, since the tail integral po- 
tentially contains highly nonlocal contributions [49-51]. 
As such, it is probably true only in a very particular set 
of situations, such as, for example, an EMRI system in 
the adiabatic limit [52], in which the geodesic motion is 
periodic and the particle executes a large number of or- 
bits before deviating noticeably from the geodesic.^ Ob- 
viously, we would like the self-force to be valid in more 
general regimes — for example, in the final moments of 
plunge in the EMRI orbit. 

A more systematic method of "patching together" reg- 
ular expansions has been devised by Hinderer and Flana- 
gan [35] . They perform a two-timescale expansion of the 
Einstein equation and the equation of motion, tailoring 
their approach to EMRI systems. This expansion cap- 
tures both the fast dynamics of orbital motion and the 
slow dynamics of the particle's inspiral and the gravi- 
tational backreaction on the background spacetime. At 



In fact, the typical derivation of the MiSaTaQuWa equation, 
which begins with a source moving on a geodesic but ends with a 
self-force, has sometimes been called an adiabatic approximation 
[32]. However, others [49] have suggested that the substitution 

h [7] — » h [7*"'] in the equation of motion is more generally 
valid. 



each value of the "slow time," one can perform a regular 
expansion, from which the self-force can be derived as dis- 
cussed above, using the actual field and the position and 
momentum of the particle as initial data — this is one rea- 
son why the force as derived in a regular expansion should 
be expressed in terms of the actual field, rather than the 
tail integral over the entire past history of a geodesic. 
The slow dynamics are assumed to be irrelevant over the 
timescale of the regular expansion, such that the slow 
time variable appears as a fixed parameter during the 
expansion. By letting the slow time evolve continuously, 
a series of regular expansions are automatically patched 
together to arrive at a self-consistent evolution. 

Another expansion has been devised by Mino [32, 34]. 
He begins with an expansion similar to the one presented 
here, but he then performs a second expansion of each 
li ' , in such a way that each term in the expansion of /i* ' 
depends only on information from the instantaneously 
tangential worldline governed by the (n — l)th-ordcr self- 
force; that is, each term in the expansion of the leading- 
order perturbation }i ' depends only on the geodesic in- 
stantancously tangential to the true worldline, h de- 
pends only on the worldline governed by the first-order 
self-force, and so on. 

The methods developed in this paper are intended to 
complement the above approaches. It is hoped that they 
will be valid in more general contexts, though more de- 
tailed studies would be required to bear out that hope. 



III. PRELIMINARY ANALYSIS II: THE 
MOTION OF AN EXTENDED BODY 

Since the singular expansion presented in the previous 
section is based on an exact point particle source, it is 
ill-behaved beyond first order. As such, we must now 
consider methods of accounting for the extension of an 
asymptotically small body. Specifically, we must consider 
how to formulate an asymptotic expansion in which a 
representative worldline for the small body is held fixed. 

Perhaps the most obvious approach is to work with a 
body of arbitrary size and then take the limit as that size 
becomes small. Such a method has been used by Harte 
[53] in deriving self- force expressions, following the ear- 
lier work of Dixon [54]. However, in this article I will 
be interested only in approaches that treat the body as 
asymptotically small from the start. The simplest means 
of doing so is to treat the body as an effective point par- 
ticle at leading order, with finite size effects introduced 
as higher-order effective fields, as done by Galley and Hu 
[23]. However, while this approach is computationally 
efficient, allowing one to perform high-order calculations 
with (relative) ease, it requires one to introduce methods 
such as dimensional regularization and mass renormaliza- 
tion in order to arrive at meaningful results. Because of 
these undesirable requirements, I will not consider such 
a method here. 
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A. Point particle limits 

In order to move from an exactly pointlike body to an 
asymptotically small one, we must consider a family of 
metrics g{e) containing a body whose mass scales as e in 
the limit e 0. (That is, m ~ eTZ.) If each member of 
the family is to contain a body of the same type, then 
the size of the body must also approach zero with s. The 
precise scaling of the size with e is determined by the type 
of body, but this precise scaling is not generally relevant.^ 
What is relevant is the "gravitational size" — the length 
scale relevant to the metric outside the body — and this 
size always scales linearly with the mass. If the body is 
compact, as is a neutron star or a black hole, then its 
gravitational size is also its actual linear size. 

Point particle limits such as this have been used to 
derive equations of motion many times in the past, in- 
cluding in derivations of geodesic motion at leading or- 
der [5, 8, 9] and in constructing post-Newtonian limits 
[3, 55, 56]. In general, deriving corrections to geodesic 
motion requires considering two types of point particle 
limits: an outer limit, in which £ ^ at fixed coordinate 
values; and an inner limit, in which £ ^ at fixed values 
of 7~ = r/£, where r is a radial coordinate centered on the 
body. In the outer limit, the body shrinks toward zero 
size as all other distances remain roughly constant; in the 
inner limit, the small body remains a constant size while 
all other distances blow up toward infinity. Thus, the 
inner limit serves to "zoom in" on a small region around 
the body. The outer limit can be expected to be valid 
in regions where r ~ 7?,, while the inner limit can be ex- 
pected to be valid in regions where f ^ TZ, though both 
of these regions can be extended into larger domains. 

These two limits can be utilized in multiple ways. For 
example, the outer limit can be used to examine the ef- 
fect of the small body on the external spacetime, while 
the inner limit can be used to study the effect of the 
external spacetime on the metric of the small body (as 
studied in Refs. [57-59], for example). What is of inter- 
est for this paper is how the two limits mesh in a buffer 
region — a region in which m <C r <C 7?.. The metric in 
this region will determine the motion. To understand 
this, note that the buffer region, at fixed time, is approx- 
imately flat, since it is simultaneously in the asymptotic 
far zone of the body (because r 3> m) and in a small 
local patch in the external spacetime (because r <^TZ). 
Thus, in the buffer region, the linear momentum of the 
small body, or some other measure of motion, can be de- 
fined. Speaking roughly, an equation for the derivative 
of this linear momentum will then provide an equation 



However, the calculations in this paper require the existence of 
a vacuum region of radius m <^ r ^ "R, around the body. If the 
body is not sufficiently compact, then this region will not exist 
and my calculation will not apply. Likewise, my calculation fails 
when a body becomes tidally disrupted. 



of motion for the body. 

We can consider two basic methods of deriving equa- 
tions of motion in the buffer region. The first method is 
that of matched asymptotic expansions, in which one ap- 
proximation is constructed using the inner limit, another 
using the outer limit, and then any undetermined func- 
tions are fixed by insisting that the two approximations 
are equal to one another in the buffer region. The second 
method foregoes an explicit calculation of an approxima- 
tion in either the inner or outer limit (or both), instead 
working entirely in the buffer region and using some lo- 
cal definition of the motion of the body. Although both 
make use of inner and outer expansions, the two meth- 
ods are logically distinct. However, both methods have 
sometimes been referred to as the method of matched 
asymptotic expansions {e.g., in Ref. [7]). 

D'Eath was the first to apply these methods to the 
problem of motion in General Relativity. He used 
matched asymptotic expansions to show that at leading 
order, a rotating black hole moves on a geodesic of the ex- 
ternal spacetime [6, 10]. Since D'Eath's pioneering work, 
these methods have been used in many contexts: to show 
that an arbitrarily structured body follows a geodesic [7] ; 
to show that a charged body follows a worldline governed 
by the Lorentz force law [60]; to derive post-Newtonian 
equations of motion [3, 55, 56, 61, 62]; to derive general 
laws of motion due to the coupling of the body's mul- 
tipoles with those of the external spacetime [59]; and 
most pertinently, to derive the gravitational self-force 
[18, 20, 24, 28, 29]. These derivations of the gravitational 
self-force will occupy the remainder of this section. 

Let us first consider the earliest such derivation, 
performed by Mino, Sasaki, and Tanaka [18], and in 
slightly different manners by Poisson [29] and Detweiler 
[28]. These derivations take the small body to be a 
Schwarzschild black hole (with the hope that more gen- 
eral bodies would obey the same equation of motion), 
such that in the inner limit the exact metric g can be 
approximated by g = 5s(r) + H{r) + 0(£^), where 
the internal background metric 5_b is the metric of the 
isolated black hole, and H{f) consists of tidal pertur- 
bations. In the outer limit, the metric is written as 
g = g + /i[7] + 0(£^), where g is an arbitrary vacuum 
metric and ft. [7] is the perturbation due to a point par- 
ticle traveling on a worldline 7. Expanding the external 
metric in normal coordinates centered on the worldline, 
expanding the internal metric for r ^ m, and insisting 
that the results of these expansions are identical, then 
determines an equation of motion for 7. 

Mino, Sasaki, and Tanaka [25] later used a similar 
method to determine an equation of motion for a small 
Kerr black hole; they followed Thorne and Hartle's [59] 
approach of defining the spin and angular momentum of 
the body as an integral over a closed spatial surface in the 
buffer region, and they derived an equation of motion by 
combining these definitions with the assumed point par- 
ticle perturbation in the external spacetime. 

Allow me to more precisely state the underlying logic 
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of these derivations, which is as follows: Suppose there 
exists a metric g such that (1) in a region TZi, g is 
well approximated by the metric of a tidally perturbed 
black hole, (2) in a region TZ2, g is well approximated by 
the metric of some vacuum spacetime as perturbed by 
a point-like source moving on a worldline 7 C TZi, and 
(3) the regions TZi and TZ2 overlap. Then the worldline 
7 is governed by the MiSaTaQuWa equation. Alterna- 
tively, a weaker formulation might be stated as follows: 
the approximate solutions to the Einstein equation given 
hy g + h[j] and gs + H, as defined above, can be com- 
bined to form a global approximate solution if and only 
if 7 is governed by the MiSaTaQuWa equation.^ 

There are several problems with this approach. First, 
it suffers from the same problem described in the previous 
section: since the point particle solution solves the lin- 
earized Einstein equation only if the point particle travels 
on a geodesic, a non-systematic gauge relaxation must be 
invoked. Second, it does not offer any way to go beyond 
first order, since it provides no means of determining the 
external perturbations (though see Refs. [21, 63] for an 
extension to second order). Third, it is a somewhat weak 
result, with many i/s in its construction. Of course, it 
does provide an improvement over the earliest point par- 
ticle derivations, since it derives the self-force from the 
consistency of the field equation and makes no question- 
able assumptions about the behavior of singular quanti- 
ties. 

More recent derivations have been performed by Fuku- 
moto, Futamase, and Itoh [24] and Gralla and Wald 
[20]. These derivations work entirely in the buffer re- 
gion, rather than using matching; they do not assume 
that the external perturbation is that of a point particle 
at leading order; and they do not restrict the small body 
to be a Schwarzschild black hole. Fukumoto et ai, fol- 
lowing the work of Futamase [3, 55, 56] and Thornc and 
Hartle [59], defined the linear momentum of the body 
as an integral in the buffer region and derived the accel- 
eration by simply differentiating this linear momentum. 
While this derivation is quite simple relative to most oth- 
ers, it contains at least one questionable aspect: it relics 
on an assumed relationship between the body's linear 
momentum, as defined in the buffer region, and the four- 
velocity of the body's worldline (justified by an analogy 
with post-Newtonian results). 

Gralla and Wald explicitly restricted themselves to a 
regular expansion, defining the acceleration of the body 
via a regular expansion of the body's worldline, roughly 
as described in the previous section. The only question- 
able aspect of their derivation is that it writes the so- 
lution to the first-order Einstein equation as an integral 



In actuality, the method of matched asymptotic expansions, as 
it has been used in derivations of the self-force, provides signif- 
icantly weaker results than this, because at first order in e the 
relationship between points in TZi and points in 7^2 is not unique, 
which leads to a non-unique result for the acceleration [43] . 



over the past history of the leading-order, geodesic world- 
line 7'"', and it expresses the force in terms of the tail 

tail 

integral h [7'°']. As discussed above, this is not obvi- 
ously justified: in order to remain consistently within the 
domain of validity of a regular expansion, the tail should 
be cut off at some finite past time and complemented 
with an integral over an initial data surface. The prac- 
tical drawback of this derivation is that it is obviously 
limited to short timcscales, as discussed above. While 
a regular expansion such as this can be used to derive 
the self-force and then incorporated into a two-timescale 
expansion, my goal here is to provide a self-consistent 
approach in which the worldline is never treated 
geodesic. 

B. Definitions of the worldline 

While the derivations described above are increasingly 
satisfactory, none of them have satisfactorily defined the 
worldline of the asymptotically small body. To see this, 
we must examine the various definitions of this worldline. 

Most of these definitions are in terms of the outer limit 
[5, 6, 8-10, 20]. At each value of e, the body can be 
surrounded by a worldtube F^ that has a radius which 
vanishes in the limit e ^ 0; the worldline of the body 
is then defined as the limit F^^o of these worldtubes, 
which defines a curve in the limiting spacetime g(e = 0). 
Kates generalized this to the case when the limit e ^ 
is singular [7]. 

These definitions are extremely problematic, because 
the worldline seems to naturally emerge only when the 
body is exactly point-like, which is the case only in the 
limiting spacetime defined by e = 0. But in this space- 
time, since s — 0, the worldline will naturally be e- 
independent; any e-dependence would automatically be 
"pure gauge," such that the self-force could be set to zero 
over the entire domain of the regular expansion. How, 
then, can the worldline be defined such that it can accu- 
rately and meaningfully reflect the motion of the body 
for e > 0? If we reject the use of small "corrections" to a 
worldline, represented by a deviation vector, how can we 
find the worldline that the deviation vector "points" to? 

One possible definition has been suggested by Futa- 
mase [3, 24, 55, 56]. Rather than defining the worldline 
as the limit of a family of worldtubes of radius ~ e, he 
defines the worldline as the curve that remains within ev- 
ery such tube as £ 0. One can easily show that if such 
a curve exists, then it is unique. However, one can just 
as easily show that, in general, such a curve exists only 
for a short time: for any two values of e, say ei and 62, 
the interiors of the worldtubes F^^ and F^^ will intersect 
only for a brief time period t <TZ. Hence, this definition 
does not improve upon the previous one. It also has the 
disadvantage that it cannot apply to small black holes. 

Yet another approach is to define the worldline implic- 
itly. Consider how this is realized in matched asymp- 
totic expansions (see Ref. [62] for the clearest example), 
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in which the worldlinc is defined roughly as follows: if 
the perturbed metric of the external spacetime near a 
worldline 7 is equivalent (up to diffeomorphism) to the 
perturbed metric of the small body, then 7 is said to be 
the worldline of the body. In practice, this means that 
the worldline is defined operationally by the fact that, in 
the buffer region, the external metric in some coordinates 
centered on 7 is identical to the internal metric in some 
mass-centered coordinates {i.e., coordinates in which the 
internal metric has no mass dipole — see Ref. [59] for a dis- 
cussion of mass-centered coordinates defined in the buffer 
region). More precisely, suppose we are given an inner 
expansion gi{e) = gs+H in some "local" coordinates 
on a manifold Mi, and an outer expansion gsi^) = g + h 
in some "global" coordinates on a manifold Me.^ For 
example, in an EMRI consisting of a small Schwarzschild 
black hole orbiting a large Kerr black hole, the coordi- 
nates might be the Boyer-Lindquist coordinates of the 
large black hole, and the coordinates might be the 
Schwarzschild coordinates of the small black hole. We 
can always transform gE into a coordinate system (e.g., 
Fermi coordinates) centered on a worldline 7 C Me, via 
a map : Me ^ Me- The worldline 7 is then defined 
to be that of the small body if there exists a unique map 
(fi : Me ^ Mi 'm a. buffer region TZ C Me such that 
{(fi o ?/;),(a;^) = and {ip o il})t.{gE) = gi, where /* de- 
notes the push-forward of / and an equal sign indicates 
equality at the lowest common order. 

So long as one restricts one's attention to approximate 
solutions of the Einstein equations, this operational def- 
inition seems to be valid. It will undoubtedly result in 
a metric that solves the Einstein equation to some spec- 
ified order in some large region of spacetime. However, 
at first glance it might seem unlikely that the resulting 
approximate solution actually approximates any "true" 
metric (that is, any exact solution to the Einstein equa- 
tion), since it is not apparent how the e-dependent world- 
line could arise from expanding an exact solution. If we 
wish the approximate metric to be an approximation to a 
true metric, rather than just an approximate solution, we 
must at least show that the approximation could plausi- 
bly be constructed from an exact solution. And this in 
turn means that we must relate our operational definition 
of the worldline to the asymptotic behavior of a family 
of exact solutions to the Einstein equation. 

In the analysis of the point particle field equation pre- 
sented in Sec. II, I argued in favor of a singular expansion 
that holds the worldline of the particle fixed. In the case 



The subscript E is derived from the title "external spacetime" 
given to the background metric g on Me- Note that two man- 
ifolds arc gencrically required. Consider, for example, the case 
of a small black hole orbiting a large black hole. The manifold 
Ail possesses a singularity at the "position" of the small black 
hole but is otherwise smooth, while the manifold A4e possesses a 
singularity at the "position" of the large black hole but possesses 
a smooth worldline where the small black hole should be. 



of a point particle, one can easily imagine performing 
such an expansion of an exact metric, since one can eas- 
ily imagine that the exact solution could be written as a 
functional of the particle's worldline. But the situation 
seems somewhat obscure in the case of an asymptotically 
small body, for which a worldline seems to appear only 
when the size of the body becomes precisely zero. How- 
ever, the matching calculation provides some insight into 
the problem. First, although the worldline is a curve in 
the external manifold Me, need not be a curve in the 
manifold M on which the exact metric g lives; in fact, if 
the small body is a black hole, then there is obviously no 
such curve. {M might be thought of as the result of cut- 
ting out a portion of Me around the worldline 7, then 
stitching part of the manifold Mi into the excised re- 
gion.) Second, note that if the outer expansion 5-1-/1 [7], 
written in the coordinates x'^, and the inner expansion 
gs + H, written in the coordinates X^, are both approx- 
imations of the same exact solution, then they induce 
coordinates on the manifold M. This means that the 
coordinate transformation ip o ip between x^ and X^ in- 
duces a transformation between the coordinates on the 
exact spacetime. And since is defined by a curve 7, this 
coordinate transformation is parametrized by that curve, 
even though 7 is not an actual worldline in M. Thus, 
the metric in the external coordinates can naturally 
be written as a functional of 7. In this sense, the "world- 
line" is simply a function that parametrizes the metric 
(c./. the discussion in Ref. [64]). 

We can also consider this from a different angle. As 
shown by Sciama et al. [65], any exact solution of the 
Einstein equation can be written in an integral formula- 
tion. Consider a bounded vacuum region Q G M with a 
boundary dVl. At any point x in the interior of Q, the 
metric will satisfy 

g"''(x)= J ^V^^G''^''\\x,x')dS''' , (39) 

where G"'^" 1,1 (x,x') is a Green's function for the opera- 
tor 

Dp.p. = ig"''g^(pg.).^V„sv^ + R^t^p,),. (40) 

(The proof of the integral identity in Ref. [65] is restricted 
to a convex normal neighbourhood of x, but for the sake 
of argument, assume that it is valid even if f2 extends be- 
yond the normal neighbourhood.) Assume that in some 
region lA around the body, a scalar field r provides a 
measure of distance from the body. The region Li need 
not include the body itself, but should have the topology 
S'^ X [ri,r2] X [ti,t2], where r2 > ri, ^2 > ti, and S"^ is a 
spatial 2-sphere around the body. Now suppose that the 
"inner" boundary of is a timelike worldtube <Z U 
of fixed radius r = 3£^e^,Q<p<l, around the 
body. The surface is parametrized by two angles 9^ 
{A~ 1,2) and by a time t. Thus, F^ is generated by a col- 
lection of timclike curves 7(5?. : 1 1— > x°'{t,3£ ,9^). Note 
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that since ^ is arbitrary within some interval, we can use 
it interchangeably with r, implying that {t,^,0^) de- 
fines a local coordinate system X'^ near the small body. 
The collection of maps ^(of ^) thus defines the coordinate 
transformation ip between these local coordinates and the 
global coordinates x^^ . 

The metric can be written as 

g^'^^y j ^Va'C^^"' u,^\n''' dtde^de^ 

+ j ^V^^G^^'^'^'dS"' . (41) 

For small values of e, the radius ^ of the tube is also 
small, so each of the curves can be expanded about 
M = 0. If 7(5?^6() is sufficiently well behaved, this ex- 
pansion is valid even if x'^{t,^, ~ 0) does not describe 
a timelike curve (or any curve) in M. However, note 
that the integrand in the above integral will generically 
diverge at ^ = 0, since the small body's contributution 
to the metric will contain terms diverging as thus, 
the integrand itself cannot be naively expanded in pow- 
ers of ^ without carefully expanding it in powers of e at 
the same time. (Alternatively, one could perform such an 
expansion and introduce some regularization method af- 
terward.) Nevertheless, in the limit of small e, the metric 
outside of the tube will naturally be expressed as a func- 
tional of a single worldline j : t i-^ x°'{t,3^ = 0). This 
curve can be made unique by demanding that the mass 
dipole of the body vanishes, up to some desired order, 
when calculated in the local coordinates {t,^,9^). 

Based on these plausibility arguments, we can reason- 
ably believe that an exact metric for a small body could 
naturally be expressed as a functional of an e-dependent 
curve that represents the motion of the body. Hence, 
we can reasonably believe that a singular expansion in 
which the metric perturbations are treated as function- 
als of a fixed, e-dependent worldline, can approximate 
an exact metric. Actually proving that the expansion in 
this paper approximates an exact solution would presum- 
ably require a monumental effort. However, as discussed 
above, the worldline of the body is uniquely defined in 
an operational sense, and the metric that depends on it 
will provide a (hopefully uniform) approximate solution 
to the Einstein equation, whether or not it provides an 
approximation to an exact solution. 

Thus, the metric in the outer limit can be taken to be 

g{x,e) = g{x) + h{x,e;j). (42) 

In this expansion, the perturbations produced by the 
body are constructed about a fixed worldline determined 
by the particular value of e at which one seeks an ap- 
proximation. Note that because the external background 
metric g does not depend on 7. it is identical to the met- 
ric g(e = 0) in the regular limit, and the manifold A^^; on 




FIG. 1: Regular limit of a family of spacetimes. The dashed 
lines indicate a limit process in which £ — > at fixed coor- 
dinate values. As e — > 0, the worldtube of the body shrinks 
to zero size, leaving a geodesic remnant curve. The remnant 
curve lies in the manifold Mo defined by e = 0. 



which it lives is identical to the manifold A^o-^" However, 
even if one writes the perturbation h as an asymptotic 
series e/i*^' [7] -|- e^/i'^' [7] -I- then the terms in the series 
will not be equal to derivatives of g(e) with respect to e 
at e = 0, since they depend on an e-dependent worldline. 
Refer to Figs. 1 and 2 for a schematic comparison be- 
tween the approximations constructed with and without 
a fixed worldline. 

Since the the inner limit is constructed in a local co- 
ordinate system around the body, which is effectively al- 
ready centered on the worldline, it need not hold the 
worldline fixed. Instead, the worldline in the outer limit 
serves to fix the location at which the inner expansion is 
constructed. That is, the fixed worldline ensures a fixed 
relationship between the inner and outer expansions. 



Note that if g is smooth, then it does not depend on 7. One 
can see this by referring to the formal solution (39). Taking 
part of the boundary to be a spatial surface S that intersects 
the timclike worldtube F and assuming that at leading order the 
interior of the worldtube is smooth, it follows that for the leading- 
order solution outside the tube, the integral over the tube can be 
replaced by an integral over a spacelike "cap" that joins smoothly 
with E. Of course, since the regular limit e — > might not exist 
at the "position" r = of the small body (for example, if the 
body is a black hole), the manifold A4o may have the "remnant" 
worldline 7(e = 0) removed from it; but this discontinuity is 
obviously removable. 
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FIG. 2: Singular limits of a family of spacetimes. The dashed 
lines correspond to the external limit, which lets the body 
shrink to zero size but keeps its motion fixed. The dotted 
lines correspond to the internal limit, which keeps the size 
of the body fixed. The worldline lies in the manifold Mo, 
but it does not correspond to the remnant curve defined by 
the regular limit; instead, it is allowed to have e-dependence, 
and it is determined by the particular value of e (in this case, 
e — e*) at which an approximate solution is sought. 



C. Outline of the construction of a uniform and 
self-consistent approximation scheme 

The foregoing discussions have made two points clear: 
First, every derivation of the gravitational self- force has 
at least one questionable aspect. Some of these ques- 
tionable aspects are fundamental — e.g., a reliance on an 
exact point particle source or an assumed form for the 
force — while others are relatively innocuous. However, 
at this point in time, nearly a dozen derivations have 
arrived at the same expression for the force — up to the 
ambiguity of whether the tail integral is to be evaluated 
over a geodesic or over the true worldline, which is not 
always clearly stated in the derivations. Thus, there can 
be little doubt, if there ever was, that the equation for 
the self-force is essentially correct. 

Second, the discussions above have clearly shown that 
the heart of the problem lies in singular perturbation 
theory, which is signalled by the failure of a regular ex- 
pansion due to the appearance of multiple length scales 
(the size ~ e of the body, the radius of curvature ^ e° of 
the external spacetimc, and the radiation-reaction time 
~ 1/e over which the worldline deviates from geodesic 
motion). This is closely related to the ambiguity men- 
tioned above: whether or not the tail integral is to be 
evaluated over the true past history of the body, or over 
a fictitious geodesic past. Straightforward analysis sug- 
gests that the integral must be over a geodesic, even if, 



contradictorily, the motion is accelerated — but that is 
a faulty conclusion based on regular perturbation the- 
ory, which is valid only for short times, and which hence 
should never have included an integral over the entire 
past. Given these facts, the most obvious way to arrive 
at a self-consistent solution is to make use of singular 
perturbation techniques. 

The derivation presented in the remainder of this paper 
is not intended to remedy the first condition mentioned 
above: it is not without questionable aspects of its own 
(though these are relatively few). What the derivation 
is intended to do is utilize singular perturbation theory 
to construct a self-consistent approximation scheme. Be- 
cause of its self-consistency, my scheme potentially pro- 
vides a uniform approximation valid over times t < TZ/s. 
Its essential feature, which distinguishes it from previous 
methods, is that it uses expansions in which the world- 
line of the object is held fixed; while this idea was taken 
for granted in some earlier derivations, it has not been 
considered explicitly before now. Within the context of 
this overall approach, I consider two specific expansions: 
an inner expansion accurate at distances r ~ £ from the 
body, and an outer expansion accurate at distances r ^ 1. 

I will now outline the structure of my approximation 
scheme. I consider a family of metrics g{s), where s > 0, 
in a large vacuum region fl outside the body. Eventually, 
the parameter e will be identified with the mass mo of 
the body at some initial time. I will ensure, by construc- 
tion, that the expansion is an asymptotic solution to the 
Einstein equation; I will hope, based on the plausibility 
arguments offered in the previous section, that the ex- 
pansion is also an asymptotic approximation to an exact 
solution. As a technicality, I assume that all quantities 
have been rescaled by the infimum of the external length 
scales in fl, such that we can meaningfully speak of the 
mass of the body, or a radial coordinate near the body, 
being small or large relative to unity. 

I choose n to lie outside a worldtube F surrounding 
the body. The tube's radius ^ is chosen to satisfy e ^ 
^ <C 1; in other words, F is chosen to be embedded in 
the buffer region where both inner and outer expansions 
are valid. Hence, from the point of view of the outer 
expansion, the radius of the worldtube is asymptotically 
small {,'% ^ 1), and its interior forms part of a smooth 
manifold A^^;, on which the external background metric 
is defined. The worldline lies in A4e, at the center of 
this smooth interior. But from the point of view of the 
inner expansion, the radius of the tube is asymptotically 
large (^ ^ e), and its interior is a subset of a manifold 
Mi, on which the internal background metric is defined, 
and in which there is potentially a black hole and no 
meaningful worldline. The worldline 7 at the center of 
the worldtube's interior — in Me — is defined to be the 
body's worldline if the body also lies at the center of the 
worldtube's interior — in A4i — in the sense that its mass 
dipole vanishes on the worldtube. Using the worldtube 
to divide the spacetime into an inner region and an outer 
region in this way serves to "cut out" the singularities 
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that would appear in the metric perturbation in the outer 
hmit, were it extended into the interior of the worldtube. 

Although I am interested in the solution outside the 
tube, I will require some information from the metric in 
the inner limit. I assume the existence of some local polar 
coordinates X" = {T, R,Q^), such that the metric can 
be expanded for e ^ while holding R = R/e, and 
T fixed. This leads to the ansatz 

g{X, e) = geiT, R, Q^) + H{T, R, Q^, e), (43) 

where H at fixed (T, R, Q^) is a perturbation beginning 
at order e. The leading-order term gsiT, R, 0^) at fixed 
T is the metric of the small body if it were isolated. For 
example, if the body is a small Schwarzschild black hole 
of ADM mass m(T). then in Schwarzschild coordinates 
5s(T,i?,e^) is given by 

ds^ = - (l- 2m{T) /R}j dT'^ + (l " 2m(r) /R^ ^ e^dR^ 
+ £^^2 (de2 + sin2ed$2) , (44) 

where m{T) = m{T)/e. Since the metric becomes one- 
dimensional at e = 0, the limit e — > is singular. As 
discussed by D'Eath [6, 10] (see also Ref. [20]), the limit 
can be made regular by rescaling time as well, such that 
T = {T — To)/s, and then rescaling the entire metric by 
a conformal factor l/e^. This is equivalent to using the 
above singular expansion and assuming that the metric 
gs and its perturbations are quasistatic (evolving only 
on timescales ~ 1). Both are equivalent to assuming that 
the exact metric contains no high-frequency oscillations 
occuring on the body's natural timescale ^ e. If one 
were interested in the effect of the external spacetime on 
the metric near the small body, one could determine the 
perturbation H by expanding it in powers of e, solving 
the perturbative Einstein equation, and then matching 
the result to the external solution. 

However, in this paper, the inner expansion will be 
used only to provide data for the outer expansion. In the 
outer limit, I expand for e ^ while holding fixed some 
global coordinates as well as the worldlinc 7. This 
leads to the ansatz 

gix,e) g{x) + h{x,e;-/), (45) 

where 

Ne 

h{x, e; 7) - E ^''^E i^; 7) + 0(e^-+^) • (46) 

n=l 

In order to solve the Einstein equation with a fixed 
worldline, I assume that the Lorenz gauge can be imposed 
everywhere in Q on the entircy of h, such that Lfj_[h] = 
0.^^ With this gauge condition, the vacuum Einstein 



Note that this is a stronger assumption than in the point particle 



equation R^^ = is reduced to a weakly nonlinear wave 
equation that can be expanded and solved at fixed 7, 
leading to the sequence of wave equations 

E,4h'E]^0, (47) 
E^,[h'E]=26^R^,[hlEl (48) 

and so on, where E^^, is the wave operator defined in 
Eq. (6). I discuss the formal solution to these equations 
in Sec. V. 

For simplicity. I assume that each term in the expan- 
sion of the metric perturbation minimally violates the 
Lorenz gauge, in the sense that if a solution truncated 
at some finite order violates the Lorenz gauge, then that 
violation is solely due to the acceleration. Again solving 
at fixed 7; this assumption leads to the equations 

ii:"[C]=0, (49) 

L'^'[h'E]--L^^'[hE], (50) 

which follow from an assumed expansion of the accelera- 
tion: 

a,(t, e) = a^\t) + ea'l\t-^) + 0{e^) . (51) 

I remind the reader that Ln is the gauge operator defined 
in Eq. (7), L^'^if] = ml^^,,,, and L<^'[/] consists of 
the terms in L[f] that arc linear in a'^'. 

In Sec. IV, following the approach of Kates [7] and 
Gralla and Wald [20], I determine the acceleration by 
solving the Einstein equation in the buffer region around 
the body. In this region, I work in Fermi normal co- 
ordinates centered on the worldlinc, and I expand the 
metric for small e and r, allowing each coefficient to have 
a functional dependence on 7. I never make use of the 
global coordinates x'^ or the local coordinates X'^, as- 
suming only that they can be transformed into Fermi co- 
ordinates in the buffer region — which is necessarily true 
if both the inner and outer expansions are approxima- 
tions to the same exact solution. Since the worldline is 
fixed, its acceleration will appear explicitly in the exter- 
nal background metric in Fermi coordinates. Hence, by 
solving the Einstein equation in these coordinates, the 
acceleration will naturally be determined. Although I 
perform this calculation in the Lorenz gauge, the choice 
of gauge should be of little significance. 

The result of this calculation is that the external metric 
perturbation in the buffer region is expressed as the sum 
of two solutions: one solution that diverges at r = and 



case, because if the metric is given in some other gauge, the gauge 
vector(s) transforming to the Lorenz gauge must satisfy not only 
some weakly nonlinear wave equation, but also some suitable 
boundary conditions on the worldtube F. However, in practice 
I will be satisfied by the existence of an approximate solution 
to the Einstein equation that approximately satisfies the gauge 
condition up to errors of order e^. 
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which is entirely determined from a combination of (i) 
the multipole moments of the internal background metric 
qb, (ii) the Riemann tensor of the external background g, 
and (iii) the acceleration of the worldline 7; and a second 
solution that is regular at r = and must be determined 
from the global past history of the body. At leading 
order, these two solutions arc identified as the Detweiler- 
Whiting singular and regular fields and [22], and 
the self- force is determined entirely by . Along with 
the self- force, the acceleration of the worldline includes 
the Papapetrou spin-force. This leaves us with the self- 
force in terms of the metric perturbation induced by the 
body. 

In Sec. V, I proceed to obtain a global, formal solution 
for the metric perturbation in the Lorcnz gauge. Fol- 
lowing the method of D'Eath [6, 10], I write the formal 
solution to the wave equation in an integral representa- 
tion, whereby the value of the metric perturbation at any 
point in the exterior region is related to an integral over 
the worldtube around the body. Since the tube is chosen 
to lie in the buffer region, the previously obtained expan- 
sion in that region then serves to provide the boundary 
data on the tube. This approach allows me to deter- 
mine in the bufi^er-region expansion by appealing to 
the consistency of the integral representation of the wave 
equation. Given the results of the buffer-region expan- 
sion as boundary values, evaluating the integral represen- 
tation at a point just outside the worldtube must return 
the general solution in the buffer region. This consistency 
condition determines the unknown functions in terms of 
a tail integral. With the solution in the buffer region de- 
termined, the worldline is also determined; at the same 
time, since the boundary values are determined, the so- 
lution in the external spacetime is also determined. 



IV. THE SELF-FORCE IN TERMS OF THE 
METRIC PERTURBATION 

A. Expansion of the metric in the buffer region 

In the buffer region, I adopt Fermi normal coordinates 
(i, centered on the worldline 7. This coordinate sys- 
tem is constructed from a tetrad e" , where on the world- 
line Cq = u" and e" is a spatial triad, and off the world- 
line the tetrad is defined by parallel propagation along a 
spatial geodesic perpendicular to u". Refer to Rcf. [29] 
for a detailed description. 

Although the solution to the wave equation is more 
naturally expressed in terms of retarded coordinates [29] , 
in the calculations here Fermi coordinates are more ad- 
vantageous; for example, the solution to the wave equa- 
tion with a point particle source is expressed as an in- 
tegral over the worldline up to a retarded time tret, but 
in my calculation, the solution to the wave equation will 
be expressed as an integral over a worldtube, which will 
be evaluated just as easily in Fermi coordinates as in 
retarded coordinates. Thus, in this calculation, the sim- 



pler form of the background metric in Fermi coordinates 
outweighs the advantages of retarded coordinates. 

I will be interested only in components in the Cartesian 
coordinates {t,x°'), but I will express these components 
in terms of the geodesic distance r = y/ SijX'^x^ and the 
angles 9^, which are defined in the usual way in terms of 
x". I also introduce the unit one-form Ua = daV, which 
has components Ua = (0, Ua) in Fermi coordinates, where 
Ua = 5abX^ /r. Note that this one-form has the convenient 
property that rf = g°'^np = (0, n"), where = 5°'^nii. I 
will use the multi-index notation = n'^ ...n** = n*^ - '* . 
Finally, I define the coordinate one-forms 1^ = dat and 
x° = dax''. 

Since r ^ 1 in the buffer region, both the background 
metric g and the external perturbations lij can be ex- 
panded for small r. In particular, the components of the 
background metric are given by the following standard 
result: 

gtt = -(1 + m,7i*)2 - r^Ro^ojn'^ + ©(r^), (52) 
3ta = -^r2i?o„,7i^J+0(r3), (53) 

gab = Sab - lr^Ra^bJn'^ + 0(^3). (54) 

Here the components of the Riemann tensor are evalu- 
ated on the worldline, and are therefore functions of t 
(and potentially 7) only. For later convenience, I define 
the electric-type tidal field 

Sab = RaObO, (55) 

and the magnetic-type tidal field 

13ab = \^a'''R^bcd- (56) 

Each of these fields is symmetric tracc-frec (STF) with 
respect to the Euclidean spatial metric 5ab- Solving the 
Einstein equation in the buffer region will determine an 
analogous expansion for the perturbations. 

However, before beginning that calculation, one should 
note that the coordinate transformation x"(f,x°) be- 
tween Fermi coordinates and the global coordinates is 
e-dependent, since Fermi coordinates are tethered to an 
e-dependent worldline. If one were using a regular ex- 
pansion, then this coordinate transformation would de- 
volve into a background coordinate transformation to a 
Fermi coordinate system centered on a geodesic world- 
line, combined with a gauge transformation to account 
for the e-dependcncc. But in the present singular expan- 
sion, the transformation is purely a background transfor- 
mation, because the e-dependence in the transformation 
is reducible to the e-dependence in the fixed worldline. 

The transformation hence induces not only new e- 
dependence into the perturbations lij^ , but also e- 
dependence in the background metric g. This new 
e-dependence takes two forms: a functional depen- 
dence on z"-{t) ~ x°'{t,x''' = 0), the coordinate form 
of the worldline written in the global coordinates x"; 
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and a dependence on the acceleration vector a"{t) on 
that worldline. For example, the first type of depen- 
dence appears in the components of the Riemann ten- 
sor in Fermi coordinates, which are related to the com- 
ponents in the global coordinates via the relationship 
RijKL{t) = Rap-ysiz'^ {t))ef Cjc'J^ej^. The second type of 
e-dependence consists of factors of the acceleration a'' (t) , 
which has the assumed expansion given in Eq. (51). 

Hence, in the buffer region we can opt to work with the 
quantities g and He, which are defined with a fixed, or 
we can opt to reexpand these quantities by substituting 
into them the expansion of a. (In either case, we would 
still hold fixed the functional dependence on z'^.) Sub- 
stituting the expansion of a in Fermi coordinates yields 
the buffer-region expansions 

g^.=gl"l{t,x'';-f) + 0{e), (57) 
= Ci(<,a;";7) + 0(e), (58) 

where g"" = g\^^^io) and ft.'"' = ft'^1a=a(o)- solving 
the Einstein equation, I will make use of both the original 

quantities g and lij , and the buffer-region quantities g'"' 

and Ii ' . 

Now, in order to determine a solution to the Einstein 
equation in the buffer region, I must first determine the 
general form of an expansion in powers of r for the met- 
ric perturbations lij. To accomplish this, I consider the 
form of the internal metric gs + H. I assume that in 
the buffer region there exists a smooth coordinate trans- 
formation between the local coordinates (T, R, 8^) and 
the Fermi coordinates (t, a;°) such that T ^ t, R ^ r, 
and 0'^ ~ 9^. The buffer region corresponds to asymp- 
totic infinity r ^ s (or f ^ 1) in the internal spacetime. 
So after re-expressing f as r/e, the internal background 
metric can be expanded as 

gBo.f^{t,f,0^) = a'sipit^o^). (59) 

n>0 

There is no a priori reason to exclude negative values of n, 
since gs is an unknown function of f. However, since the 
internal and external solutions must be approximations 
to the same metric, they must agree with one another. 
And since the external expansion has no negative powers 
of e, neither has the internal expansion. Furthermore, 
since g + h = gs + H, we must have g^' = gi^"^ = 
0), since these are the only terms independent of both 
£ and r. Thus, noting that g(a;" = 0) = 77, where rj = 
diag(— 1, 1,1,1), I can write 

+ ^)\%Ui^0^) + O{eVr^), (60) 

implying that the internal background spacetime is 
asymptotically flat. 



I assume that the perturbation H can be similarly ex- 
panded in powers of e at fixed f, 

+ e^H2l{t,f,e^-^) + 0{e^), (61) 

and that each coefficient can be expanded in powers of 
1/f = e/r to yield 




+ 0{e^lr^) , (62) 

+ eHnrH^:^'^'"^+0{syr), (63) 
e^H^:iif) = 0{s^e\er^,r'), (64) 

where iJ<""'', the coefficient of e" and r™, is a function of 
t and 9^ (and potentially a functional of 7). Again, the 
form of this expansion is constrained by the fact that no 
negative powers of e can appear in the buffer region. 
Note that explicit powers of r appear because ef = r. 
Also note that I allow for a logarithmic term at second 
order in e; this term arises because the retarded time in 
the internal background includes a logarithmic correction 
of the form elnr (e.g., t — r ^ i — ?'* in Schwarzschild 
coordinates). Since I seek solutions to a wave equation, 
this correction to the characteristic curves induces a cor- 
responding correction to the first-order perturbations. 

The expansion of H may or may not hold the acceler- 
ation fixed. Regardless of this choice, the general form 
of the expansion remains valid: incorporating the expan- 
sion of the acceleration would merely shuffle terms from 
one coefficient to another. And since the internal met- 
ric gs + H must equal the external metric g + h, the 
general form of the above expansions of the gs and H 
completely determines the general form of the external 
perturbations: 

h-Eap = -n-Eap + n-Eaf) + ^h^p + O [r ) , (65) 

(2) ^ 1 , (2.-2) 1 , (2,-1) , (2.0) (2,0.1.1) 

"■_Ea/3 = '^'>'EaP + ~'^Eal3 + "_Eq/3 + ^'^^'''EaP 

+ 0(r) , (66) 

where each li^ ' depends only on t and 0^, along with 
an implicit functional dependence on 7. If the internal 
expansion is performed with a held fixed, then the in- 
ternal and external quantities are related order- by-order: 

c-g-, E™^f'"-" = 5, ^T\= and = 

Since I am not concerned with determining the internal 



One might think that terms with negative powers of e could be 
allowed in the expansion of gg if they are exactly cancelled by 
terms in the expansion of H, but the differing powers of r in the 
two expansion makes this impossible. 
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perturbations, the only such relationship of interest is 
li^ ' = gg\ This equality tells us that the most di- 
vergent, r~" piece of the rtth-order perturbation Ji^ is 
defined entirely by the nth-order piece of the internal 
background metric gs, which is the metric of the body if 
it were isolated. 

To obtain a general solution to the Einstein equation, 
I write each fi^ ' as an expansion in terms of irreducible 
symmetric trace-free pieces: 

C =E<'"'^^ (67) 
e>o 

+ E Kr^if^'-' + ^a.^^iT^n''^-'] . (68) 
f.>l 

f^Eab = Sab E ^^l"-"'""^ + E '^L ""'"ab"^ 

i>o e>o 

£>1 

(69) 

Here a hat indicates that a tensor is STF with respect 
to (Sab, angular brackets () indicate the STF combina- 
tion of enclosed indices, parentheses indicate the sym- 
metric combination of enclosed indices, and all the up- 
percase script symbols are functions of time (and poten- 
tially functionals of 7) and are STF in all their indices. 
Each term in this expansion is linearly independent of 
all the other terms. All the quantities on the right-hand 
side are flat-space Cartesian tensors; their indices can be 
raised or lowered with bah- Refer to Appendix B for more 
details about this expansion. 

Now, despite its e-dependence, g is the background 
metric of the outer expansion, and I will use it to raise 
and lower indices on h. And since the wave equations (47) 
and (48) are covariant, they must still hold in the new 
coordinate system, despite the additional e-dependence. 
Thus, both equations could be solved for arbitrary accel- 
eration in the buffer region. However, due to the length 
of the calculations involved, I will instead solve the equa- 
tions 

i?„^[4']=0, (70) 
£;l-[/.'^'] = 2<5X°M'^''']+0(e), (71) 

where i?'«>[/] = i?[/]L^^(o, and 8^R'^^\!\ = 
(5^i?[/] |^_^(o, .^^ The first equation is identical to 



In analogy with the notation used for L'"' , B'^' [/] and (5^_R*^' [/] 



Eq. (47). The second equation follows directly from sub- 
stituting Eqs. (57) and (58) into Eq. (48); in the buffer 

(2) 

region, it captures the dominant behavior of , repre- 

(2) 

sented by the approximation h , but it docs not cap- 
ture its full dependence on acceleration. If one desired 
a global second-order solution, one would solve Eq. (48), 
but for my purpose, which is to determine the first-order 
acceleration a'^', Eq. (71) will suffice. 

Unlike the wave equations, the gauge conditions (49) 
and (50) already incorporate the expansion of the accel- 
eration. As such, they are unmodified by the replacement 
of the second-order wave equation (48) with its approxi- 
mation (71). So we can write 

i^l\C\=^. (72) 

^J.^'[C]=-ir^'1' (73) 

where the first equation is identical to Eq. (49), and the 
second to Eq. (50). (The second identity holds because 
L<,«)[/i<"] = L(,°'[/i*^'], since /i*" diflFers from /i*^' by a^^' 
and higher acceleration terms, which are set to zero in 
L'"'.) I remind the reader that while this gauge choice 
is important for finding the external perturbations glob- 
ally, any other choice would suffice in the buffer region 
calculation. For example, one could expand g and h in 
buffer region expansions that incorporate the expansion 
of the acceleration, enabling one to solve the full Ein- 
stein equation order- by-order in e; it might be difficult 
to make this mesh with a global expansion in the external 
spacetime, but it would suffice to determine the acceler- 
ation. Alternatively, one could construct a two-timcscale 
expansion in the buffer region, which would mesh with a 
global two-timcscale expansion of the Einstein equation 
in the external spacetime. 

As a final, important point, I assume the partial time- 
derivative of any term in an expansion is of the same 
order as the term itself. In what follows, the reader 
may safely assume that all calculations arc lengthy unless 
noted otherwise. 



B. First-order solution in the buffer region 

In principle, solving the first-order Einstein equation 
in the buffer region is straightforward. One need simply 
substitute the expansion of ft,y, given in Eq. (65), into 
the linearized wave equation (70) and the gauge condi- 
tion (72). Equating powers of r in the resulting expan- 
sions then yields a sequence of equations that can be 
solved for successively higher-order terms in lij^ . Solv- 
ing these equations consists primarily of expressing each 



would be linear in a'^', and (5^_R'^'[/] would be linear 

in a'^' and quadratic in a*^', and so on. For a function / ~ 1, 
L(">[/], and g^R'-'^'^\f\ correspond to the coefficients of 

in expansions in powers of e. 
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quantity in its irreducible STF form, using the decom- 
positions (B3) and (B7); since the terms in this STF de- 
composition are hnearly independent, we can solve each 
equation term-by-term. This calculation is aided by the 
fact that Vj = .tJ^9q -t- 0(r°), so for example, the wave 
operator Eap consists of a flat-space Laplacian d°'da plus 
corrections of order l/r. Appendix B also contains many 
useful identities, particularly dar = Ua, n^dafi^ = 0, 
and the fact that is an eigenvector of the flat-space 
Laplacian: i.e., d°'dan^ ~ ~^^^^r^n^ ■ Because the cal- 
culation consists mostly of simple, albeit lengthy algebra, 
I will for the most part simply summarize results. 

Of course, the Einstein equation in the buffer region 
does not completely determine the solution: auxiliary 
boundary data must also be provided. Since the most 
singular term, fi^^j^, is the ordcr-l/r term in the inter- 
nal background metric gs, it will be fully determined in 
terms of the mass of the internal spacetime. Some of the 
subleading terms will also be determined by the mass, 
while others will remain unknown. The unknowns form 
the Detweiler- Whiting regular field; they will eventually 
be expressed in terms of a tail integral in Sec. V. 

So, we begin with the the most divergent term in the 
wave equation: the order-l/r^, flat-space Laplacian term 



-d'dJ'^Z = 0- (74) 



The tt-component of this equation is 



= -^£(^+l)^^''-''n^, (75) 

£>0 



from which we read off that is arbitrary and 

must vanish for all £ > 1. The to-component is 



= -^(£ + l)(^ + 2)^<^^'-^'n/ 

-^^(^-l)^^->n^-i 

e>i 

l)eabc^','L-if^b'-\ (76) 

£>1 



from which we read off that ^a^' ^^' is arbitrary and all 
other coefficients must vanish. Lastly, the afo-component 



is 

-^(^ + 2)(£ + 3)4^'-^'n,fc^ 

-^^(£-fl)^-i;,n,^^-i 

i>i 

-^(^+l)(£ + 2)e,,(,ri,)^^-i?#<^-:\ 
e>i 

_^(^_2)(£-l)^J-^,n^-2 

e>2 

from which we read off that J^'''^'"^' and are ar- 

bitrary and all other coefficients must vanish. Thus, we 
find that the wave equation constrains Ii^ ' to be 

+ (<5afc^<^'-^'+=^X"')<4- (78) 

This is further constrained by the most divergent, 1/r^ 
term in the gauge condition: 

-^^eIW + ^"a?7^"Ci = 0- (79) 

From the ^-component of this equation, we read off 
= 0; from the a-component, JT'^ = ^Z'^-^' and 
^ 0. Thus, 

^Ea0 depends only on a single func- 
tion of time, ^/'■'-•''-K By the definition of ADM mass, this 
function (times e) must be twice the mass of the internal 

background spacetime. Thus, fi^ ' is fully determined 
to be 

f^sZ = '^Mt){tc.tp + 5a6<4), (80) 

where m(t) is defined to be the mass at time t divided 
by the initial mass e = mo- (Alternatively, we could set 
£ equal to unity at the end of the calculation, in which 
case m would simply be the mass at time t; obviously, the 
difference between the two approaches is immaterial.) 

At the next order, fi^ ' , along with the acceleration of 
the worldline and the time-derivative of the mass, first 
appears in the Einstein equation. The order-l/r^ term 
in the wave equation is 

d^dc^sZ = -'^acn''{3tatp - SabX^Xp), (81) 

where the terms on the right arise from □ acting on 
^Ii^ This equation constrains ' to be 

CI = ^r', (82) 
f^Eai = Sab (^<^'"' - ma^n^) + Jf^^'^K 
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Substituting this result into the order- l/r term in the 
gauge condition, we find 



4 Am 

—ta,dtm+ a^'"^ = 0. 

r r 



(83) 



Thus, both the leading-order part of the acceleration and 
the rate of change of the mass of the body vanish: 



— = a<°' = 
dt ' * 



(84) 



At the next order, rh^ , along with squares and 
derivatives of the acceleration, first appear in the Ein- 
stein equation, and the tidal fields of the external back- 
ground couple to -ti^ ^' . The order- 1/ r term in the wave 
equation becomes 



oro 2\ ,(1.1) 20m „ 3m 
ro OcA — I hEtt — 5 — '^y" a/^iaj-^n •' 



3r 



8m , 

H a,;a , 

r 



(85) 



2\ ,(1,1) 8m 
r . 



rd^d, + 3 ) = ----ea^.Bife'' - (86) 



4m , 
r 



2\ (1,1) 20m, 76m 

rd dc + -] h^ab = -^SabCijU ^^ah 

r J 9r 9r 



16m j ^ 8m 
—t^^nb)^ + —a^aab) 



+ ^(5a6(|aia'- 3a(,aj)n*^) . (87) 



From the tt-component, we read off that is arbi- 



3)^ 



trary, i/*^'^' = Ama^a', and = | 

from the ta-component, "^/j^'^', and ^j'^'^' are ar- 

bitrary, = — 2mdi, and ^-j^'^' = ^mBij] from 

the a6 component, J(f^^^'^\ '^'^ijk^ ^-iid ^^ii 

arbitrary, and = Ima^a', ^j-^'^' = ~^m£ij + 

|ma(iaj), ^Z^'''' = fmfy, and ^j'''' = --^mfij + 

Substituting this into the order-r" terms in the gauge 
condition, we find 

0=in^ + ra')C'. " h'^'ina - rda)}iE'xl 



- OthEat - 2^/ OthEf^^ta 



where the equation is to be evaluated at a = a'"' = 0. 
From the ^-component, we read off 

^<i'i' = idt (i/*''"' + 3jr<'-''>) . (89) 

From the a-componcnt, 

,^(1.1) _ 3 / .vy(i.i) ^(1,1) 



a 10 



It is understood that both these equations hold only when 
evaluated at a = a*"* . 



Thus, the order-r component of is 



(1.1) 

^Ett 



,(1,1) 



: Amaio} -f -I- |m£i, 

-I- |ma(iaj)n'^ 



(91) 



_ ^'^'^'n, - 2mda + "^i-'^'n* + ea^'&f■''n' 
^slb = 5ab{\'maia} + J^^'' ''n* - ^mSijff^ 

where and ,j^a'^^ ^'^^ constrained to satisfy 

Eqs. (89) and (90). 

To summarize the results of this section, we have 



^Ec^p = 7^Eai3 + + ^^Eo^p + 0(r ), where h^^^ 

is given in Eq. (80), h^^g^ is given in Eq. (82), and h^^p 
is given in Eq. (91). In addition, we have determined 
that the ADM mass of the internal background space- 
time is time- independent, and that the acceleration of 
the body's worldline vanishes at leading order. 



C. Second-order solution in the buffer region 

Though the calculations are much lengthier, solving 
the second-order Einstein equation in the buffer region 
is essentially no different than solving the first. I seek to 
solve the approximate wave equation (71), along with the 
gauge condition (73), for the second-order perturbation 

determine the acceleration a'^'. 

,(0) _ 



doing so will also, more importantly. 
In this calculation, the 
acceleration is set to a = a'"' = everywhere except in 
the left-hand side of the gauge condition, L<^'[ft,y], which 
is linear in a'^*. 

Substituting the expansion 



(2) 1 (2,-2) 1 (2,-1) (2,0) (2,0, In) 

'^ap = :pi^ap +-Kp +K!3 +Mr)Ki3 

+ 0{e,r) 



(92) 



and the results for from the previous section into the 
wave equation and the gauge condition again yields a 
sequence of equations that can be solved for coefficients 
of successively higher-order powers (and logarithms) of 
r. Due to its length, the expansion of the second-order 
Ricci tensor is given in Appendix D. Note that since 
the approximate wave equation (71) contains an explicit 

(2) 

0{e) correction, h will be determined only up to 0{e) 
corrections. For simplicity, I omit these 0(e) symbols 
from the equations in this section; note, however, that 
these corrections do not effect the gauge condition, as 
discussed above. 
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To begin, the most divergent, order-l/r^ term in the 
wave equation reads 



1 



Am? , , 

- —to^tp, (93) 



where the right-hand side is the most divergent part 
of the second-order Ricci tensor, as given in Eq. (D3). 
From the tt-component of this equation, we read off 
_g^(2,-2) _ _277^2^ i^jjg^i; s/^^'-^'' is arbitrary. From the 
ta-component, and are arbitrary. 

From the afo-component, = |m^, if*^"^' = —7m?, 

and ^a''' ' ^abc "' and J^^'^'' are arbitrary. 
The most divergent, order- 1/ terms in the gauge con- 

(2,-2) 

dition similarly involve only h ; they read 



ADM angular momentum. Mi is a mass dipole term; it 
is what would result from a transformation x"" x"' + 
M°' /m applied to the 1/r term in . 5^ is a spin dipole 
term. Thus, the order-l/r^ part of li ' reads 

" = -2m^ + 2M,.n 



,(2,-2) 
\t 

,(2,-2) 

'''ta 

,(2,-2) 



(98) 



At the next order, 1/r'^, because the acceleration is set 
to zero, ii ' does not contribute to E^°^[}i '], and li ' 
does not contribute to 6^R'°''[h ]. The wave equation 
hence reads 



1 .<^'-i' 



A2 n(0.-3) 



(99) 



1 (r9^ - 2^") /I'^f ' - ^^'"^ (^^« - 2n,)/i';:'' = 0. 

(94) 

After substituting the results from the wave equation, the 
t-component of this equation determines that = 0. 

The a-component determines that -^abc^^ ~ ^' ■^ab ^^ ~ 
0, and 



Thus, the order-l/r^ part of h is given by 



(95) 



where 5''R'°f' 



is given in Eqs. (D4)-(D6). The 

2QCJ 



tt-component of this equation implies r^d'^dc^itt ' - 
GmJ^'^'^'fi*-', from which we read off that jz/'^ "^' is ar- 



The ta-component im- 



bitrary and ■s^ij'^^ 

plies r^d'^dcfita ' = 6m'^/^'"'n^, from which we read off 
-1) = -m^/'-"' and that '^^^•-'^ is arbitrary. The ab- 



9(1,0) 



component implies 



r^d^drK 



(2,-1) 

c'l'ab 



6m (j^< 



+ 2m5aby4.""''f^"\ 



(100) 



- 7 flab + ^^'^"""nj). 



(96) 



from which we read off that is arbitrary, J^^^^'"^' 

i(i^'^-"> + .-^(^-'»), J?f^-'' 

2mJ^'j^'°', and ■J^^b'^^ is arbitrary. This restricts /i^ 
to the form 



9(1.0) a<^(2,-l) 



where is given by Eq. (95). 

The metric perturbation in this form depends on five 
free functions of time. However, from calculations in flat 
spacetime, we know that order-e^/r^ terms in the metric 
perturbation can be written in terms of two free func- 
tions: a mass dipole and a spin dipole. We transform the 
perturbation into this "canonical" form by performing 
a gauge transformation (c./. Ref. [66]). The transfor- 
mation is generated by = —^^''^''^-'ta — ^■^a'^'^^^a^ 
the effect of which is to remove and '"^a'^^ from 

the metric. This transformation is a refinement of the 
Lorenz gauge. (Effects at higher order in e and r will be 
automatically incorporated into the higher-order pertur- 
bations.) The condition - 3J^;<"'-'' + 3i/f = 
then becomes J^a'^''^^ = -c/J^-"^' . The remaining two fmic- 
tions are related to the ADM momenta of the internal 
spacetime: 



2M, 



2S„ 



(97) 



where Mi is such that dtMi is proportional to the ADM 
linear momentum of the internal spacetime, and Si is the 



C = Sab - im^<.^'«'n'^) (101) 

-m(^<^''" -h^<^-"')nafc 

(2-2) (2-1) 

We next substitute /i ' and /i into the order- l/r'^ 
terms in the gauge condition. The t-component becomes 



1 



(4m'? 



+ UdtM, + 3<^f = 0, 



from which we read off 



-4aM, - Im-^' 



(1,0) 



(102) 
(103) 



And the a-component becomes 
1 



= 



^(1,0) _ 1_^(2,-1) ^ 1^(2,-1)) 

- 2e„,n^5t5^ (104) 
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from which wc read off 



(106) 



and that the angular momentum of the internal back- 
ground is constant at leading order: 



dtS' = 0. 



(107) 



(2) 

Thus, the order-l/r term in h is given by 



^m(^'^'"> + jr<^'"') 



(108) 



-m(^<^'«'-|-jr'^''")"a6 

+ 2mJ^^->n,/ + |m.y^J,^-'. 



Note a peculiar feature of this term: the undetermined 
function appears in precisely the form of a mass 

monopole. The value of this function will never be deter- 
mined (though its time-dependence will be). This am- 
biguity arises because the mass m that I have defined is 
the mass of the internal background spacetime, which is 
based on the internal limit process that holds e/r fixed. 
A term of the form e" /r appears as a perturbation of this 
background, even when, as in this case, it is part of the 
mass monopole of the body. This is equivalent to the 
ambiguity in any expansion in one's choice of small pa- 
rameter: one could expand in powers of £, or one could 
expand in powers of £-|-£^, and so on. It is also equivalent 
to the ambiguity in defining the mass of a non-isolated 
body; whether the "mass" of the body is taken to be m 
or m -I- is a matter of taste. As we shall dis- 

cover, the time-dependent part of J^^'^ ^^' is constructed 
from the tail terms in the first-order metric perturbation. 
Hence, the ambiguity in the definition of the mass is, at 
least in part, equivalent to whether or not one chooses to 
include the free gravitational field induced by the body 
in what one calls its mass. In any case, I will define the 
"correction" to the mass as 5m = 

We next move to the order-ln(r) /r^ terms in the wave 
equation, and the order-ln(r)/r terms in the gauge con- 
dition, which read 



Inr (2,0, i„) 
—d dch^f} = 0, 



2'/ ■^a^a"'ab 



(109) 
(110) 



From this we determine 



(Ill) 



Finally, we arrive at the order-l/r^ terms in the wave 
equation. At this order, the body's tidal moments be- 
come coupled to those of the external background. The 
equation reads 

2 



,(2,0) 



(112) 



where Eap comprises the contributions from h and 

/i''' given in Eqs. (DIG), (D15), and (D21). The con- 
tribution from the second-order Ricci tensor is given in 
Eqs. (D7)-(D9). 

Foregoing the details, after some algebra we can read 
off the solution 

C = ^<^-°> -f ^l^-'^n' + <f "'n'^' + <^"'n^-'''= (113) 



be 



n 



^3 

2 

S»(2,0) 



rib 



(114) 



Sab (^'^-"^ + + ^jr'ri'^') 



^(2,0) 



abi 



+^'\ay^;fn,, (115) 



where each of the STF tensors is listed in Table I. 

In solving Eq. (112), we also find that the logarithmic 
term in the expansion becomes uniquely determined: 



(2,0, In) 

a/3 



m SabX'^Xp. 



(116) 



This term arises because the sources in the wave equation 
(112) contain a term (x Sab, which cannot be equated to 
any term in d'^dcfig^^f, ' . Thus, the wave equation cannot 
be satisfied without including a logarithmic term. Recall 
that the logarithmic term arises at the order we would 
expect it to: the first-order perturbation alters the null 
cone of the spacetime, such that, e.g., t — r t — r — 
2em Inr, which naturally introduces a correction ~ In r 
to the order-e terms in the solution to the wave equation. 

We now move to the final equation in the buffer re- 
gion: the order-l/r gauge condition. This condition will 

determine the acceleration a<^' 
contributes to Eq. (73): 

[4)] ^ i!!^aa)^«. (117) 

(2) 

The contribution from h is most easily calculated by 
making use of Eqs. (B24) and (B25). After some alge- 
bra, we find that the i-component of the gauge condition 
reduces to 



At this order, first 







—dtSm H dt^ 



/a-.o) 



r 3r 
and the a-component reduces to 



+ ^dtJf'''"\ (118) 



3r 







4 „9 , ^ 4to ,,, 
r r 



4 

-£a^^P 
r 



4 

-Ba,S' 
r 



^^(1,1) 



r 



(119) 
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TABLE I: Symmetric trace-free tensors appearing in the order-r'^ part of the metric perturbation h in the buffer region 
around the body. Each tensor is a function of the proper time t on the worldline 7, and each is STF with respect to the 
flat-space metric Sij. 
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The reader is reminded that these equations are vahd tion 



0, except in the 
In the following 



only when evaluated at a{t) = a*°'(i) 
term ^a^^' that arose from L'-^^ Ii^ 

subsection, this will allow me to swap partial derivatives 
with covariant derivatives on the worldline. 

The i-component determines the rate of change of the 
mass correction Sm. It can be immediately integrated to 
find 



Sm{t) = Sm{Q) + gw 



<^''"(i) + 5jr<^-°>(t)] 
im p^'^-"' (0) + S^*^'"' (0)] . (120) 



If one felt so inclined, one could incorporate (5to(0) into 
the leading-order mass m. The time-dependent terms 
correspond to the effective mass created by the gravita- 
tional waves emitted by the body. 

The a-component of the gauge condition determines 
the acceleration of the worldline. Note the most impor- 
tant feature of Eq. (119), which is that it contains two 
types of accelerations: dfMi and a*^'. The first type 
is the second time derivative of the body's mass dipole 
(or the first derivative of its ADM linear momentum), as 
measured in a frame centered on the worldline 7. The 
second type of acceleration is the covariant acceleration 
of the worldline relative to the external spacetime. In 



1^(1.1) 



(121) 



This equation of motion contains two types of terms: a 
Papapetrou spin force, given by — S'jB* , which arises due 
to the coupling of the body's spin to the local magnetic- 
type tidal field of the external spacetime; and a self-force, 
arising from homogenous terms in the wave equation. 

Note that if we had followed the path of Gralla and 
Wald [20], wc would have identified a'"' as the accelera- 
tion of the worldline 7*°'. This would be the only actual 
worldline in play; all the corrections to the motion would 
be vectors defined on it. Hence, when we found a*"' = 0, 
we would have identified the worldline as a geodesic, and 
there would be no corrections a*"' for n > 0. We would 
then have arrived at the equation of motion 



5(1,0) 



^SiB'. (122) 



This is precisely the equation of motion derived by Gralla 
and Wald. It describes the drift of the body away from 
the reference geodesic If the external background is 
flat, then the mass dipole has a valid meaning as a dis- 
placement vector regardless of its magnitude; the second 
derivative dfMi then provides a perfectly valid defini- 
tion of the body's acceleration for all times. However, if 
the external background is curved, then Mi has meaning 



other words, d^Mi corresponds to the acceleration of only if the body is "close" to the worldline. Thus, d^Mi 



the body's center of mass relative to the center of the 
coordinate system, while measures the acceleration 
of the coordinate system itself. We define the world- 
line to be that of the body if the mass dipole vanishes 
for all times, meaning that the body is centered on the 
worldline for all times. If we start with initial conditions 
Mi(0) = = dtMi{Q)^ then the mass dipole remains zero 
for all times if and only if the worldline satisfies the equa- 



ls a meaningful acceleration only for short times, since it 
will generically grow large as the body drifts away from 
the reference worldline. On that short timescale of valid- 
ity, the deviation vector defined by accurately points 
from 7'"^ to a "corrected" worldline 7; that worldhne, 
the approximate equation of motion of which is given in 
Eq. (121), accurately tracks the motion of the body. Af- 
ter a short time, when the mass dipole grows large and 



24 



the regular expansion scheme begins to break down, the 
deviation vector will no longer correctly point to the cor- 
rected worldlinc. 

To summarize the results of this section, the second- 
order perturbation in the buffer region is given by /i^^ = 

77^a/3 +f^a/3 + /^a/s + ln(r)/i„^ 0(e, r) , where 
is given in Eq. (98), li" ^' in Eq. (108), ft.'^'"' in 
Eq. (113), and in Eq. (116). At order e^jr^, the 

metric is written in terms of the mass and spin dipoles 
of the internal background metric gs- The mass dipole 
is set to zero by an appropriate choice of worldline. At 
leading order in e, the body's spin is constant along the 
worldline. At order e^/r, there arises an effective cor- 
rection to the body's mass, given by Eq. (120), and the 
order-e term in the expansion of the body's acceleration 
is given by Eq. (121). 



D. Notes on the force and the field in the buff'er 
region 

The foregoing calculation completes the derivation of 
the gravitational self-force, in the sense that, given the 
metric perturbation in the neighbourhood of the body, 
the self-force is uniquely determined by irreducible pieces 
of that perturbation. Explicitly, the terms that appear 
in the self-force are given by 

= 1^/ "aCdf^, (123) 

^i^-"' = C- (124) 

Making use of the fact that li^^^ ' is a monopole, we can 
write the acceleration as 

This is all that is needed to incorporate the motion of the 
body into a dynamical system that can be numerically 
evolved; at each timestep, one simply needs to calculate 
the field near the worldline and decompose it into irre- 
ducible pieces in order to determine the acceleration of 
the body. (Obviously, such a procedure is vastly more 
complicated than what I have just implied [67-73].) The 
remaining difficulty is to actually determine the field at 
each timestep. In the next section, I will write down 
formal expressions for the metric perturbation, and in 
particular, I will determine the metric perturbation at 
the location of the body in terms of a tail integral. 

However, before doing so, I will emphasize some im- 
portant features of the self-force and the field near the 
body. First, note that the first-order external field fij 
separates into two distinct pieces. There is the singular 



piece given by 

+ r2 (|a(,aj) + ffy) n'^ } -f 0(r2) (126) 
/if, = ~2mraa + Imrea.jBih''' + 0{r^) (127) 
h^b = ^{'5ifc[l - ^raX + Ir'^ay 

+ r'^ (|a(iaj) - JE^ij) n"^^] + 2a(aa6) 

- fr^Sab + V^lfib)?^ + 0{r^) . (128) 

This field is a solution to the homogenous wave equation 
for r > 0, but it is divergent at r = 0. It is the generaliza- 
tion of the 1 /r Newtonian field of the body, as perturbed 
by the tidal fields of the external spacetime g. Follow- 
ing the method used in Sec. 5.3.5 of Ref. [29], one can 
easily show that this is precisely the Detweiler- Whiting 
singular field, given by 

hip ^4m f Glp^,p,u'''u^'dt', (129) 

where Gf^^^,^, is the singular Green's function (defined 
in Appendix A). 

Next, there is the regular field h'^ = fij — , given 
by 

+ 0(r2), (131) 

This field is a solution to the homogeneous wave equa- 
tion even at r = 0. It is a free radiation field in the 
neighbourhood of the body. And it contains all the free 
functions in the buffer-region expansion. 

Now, the acceleration of the body is given by 

a^^^ = \dahf,-dthl-^S,Bl, (133) 

which we can rewrite as 

+ ^R'^p^&u^S^' (134) 

where S'^^ = eje^^^e'^'^^ S j . In other words, a non-spinning 
body (for which S'^^ = 0), moves on a geodesic of a space- 
time g + eh^, where is a free radiation field in the 
neighbourhood of the body; a local observer would mea- 
sure the "background spacetime," in which the body is 
in free fall, to be g + eh^, rather than g. If we performed 
a transformation into Fermi coordinates in g + eh^ , the 
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FIG. 3: The spacetime region fl is bounded by the union of 
the spacelike surface E, the timehke worldtube F, and the null 
surface J. 

metric would contain no acceleration term, and it would 
take the simple form of a smooth background plus a sin- 
gular perturbation. These points were first realized by 
Detweiler and Whiting [22] and since emphasized espe- 
cially by Detweiler [28]. They are, perhaps, made espe- 
cially clear in the derivation presented here, which natu- 
rally demarcates the singular and regular fields. 

As a final note, I remind the reader that the equation 
of motion would contain an antidamping term [19, 47, 48] 
if I had not assumed that the acceleration possesses an 
expansion of the form given in Eq. (51). 

V. THE METRIC PERTURBATION 

A. Integral formulation of the Einstein equation in 
the external spacetime 

A solution to the self-force problem consists of a pair 
(7,/i). In the previous section, we have determined the 
equation of motion of 7; we now require a means of de- 
termining the metric perturbation. Specifically, on the 
external manifold Me, I seek an approximate solution 
of Einstein's equation in a vacuum region H = 51 U 
where is a bounded, open subset oi Me- I now specify 
this region to be the future range of dependence of a sur- 
face formed by the union of a worldtube T and a spatial 
surface S. This implies that the future boundary of f2 
is a null surface J7. Refer to Fig. 3 for an illustration. 
The boundary of the domain is hence dil = F U 1/ U E. 
The worldtube F is defined by a constant Fermi radial 
coordinate distance r = M from the worldline j G Me- 
Since the tube is an artificial division of spacetime, and it 
may be located anywhere in the buffer region, any valid 
solution cannot depend on The spatial surface S is 
chosen to intersect F at the initial time t = 0. 

In f2, the Lorenz gauge is imposed on the entire pertur- 
bation h, splitting the Einstein equation into the weakly 
nonlinear wave equation 

E^p [h] = 2S^R^p [h] + 0(e3) (135) 
and the gauge condition L^j [h\ — 0. Note that if a solu- 



tion to the wave equation satisfies the gauge condition on 
dfl, then the wave equation ensures that the gauge con- 
dition is satisfied everywhere. And since I have already 
determined the equation of motion using the expansion 
in the buffer region, I will hence not be interested in the 
gauge condition here. 

As discussed in Sees. II C and IIIC, I assume the ex- 
pansion haf3{x, e; 7) = ^"'^Eapi^i l) ^^'^ arrive at the 
sequence of wave equations 

E^0[h'E]=O, (136) 
E^fs[h'E]='25^Ro.p[h'E]- (137) 

Following D'Eath [6, 10, 29], I rewrite the wave equa- 
tions as integro-differcntial equations by calculating 
E[G'"^^]h - ElhlG""^^ (where G""^^ represents the ad- 
vanced Green's function for Ef^^), integrating both sides 
of the resulting equation, making use of Stokes' law, and 
finally simplifying the result using the reciprocality rela- 
tion G^Jj,^ij{x' ,x) ~ Gaf}a' i3'{x,x'). Thc resulting equa- 
tions arc 

(138) 

,(2) 1 //^ ■y'S'u^^'' u^^^ ^ -y'S' A 7cm' 

on 

-^J G^/''S^Rys'[h'E]dV'. (139) 
n 



Alternatively, we might rewrite Eq. (135) directly: 




dn 



-j- I G^/''d^Rys' [h]dV' + 0{e^). (140) 

Note that any solution to Eq. (135) in will also satisfy 
this integro-differential equation; however, because of the 
e-dependence of the true worldline 7, not every solution 
to Eq. (135) will admit an expansion satisfying thc two 
equations (136) and (137) (though a solution to thc latter 
is obviously a solution to the former). In that sense, 
Eq. (140) is more robust than Eqs. (138) and (139). 

In any case, these integral representations all have sev- 
eral important properties in common. First, the integral 
over the boundary is, in each homogeneous so- 

lution to the wave equation, while the integral over the 
interior is an inhomogeneous solution. Second, the in- 
tegral over the boundary can be split into an integral over 



The integral over the interior will also contain homogeneous solu- 
tions. However, these will be .!^-dcpcndent, and they will exactly 
cancel corresponding .^-dependent terms in the boundary inte- 
gral. 
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the worldtube F and the spatial surface E; the contribu- 
tion of the null surface J vanishes by construction. Also 
note that x must lie in the interior of 17; an alternative 
expression must be derived if x lies on the boundary [74] . 

Furthermore, the integral representations avoid any 
divergence in the second-order solution. Comparing 
Eq. (140) to the analogous expression for a point particle, 
given in Eqs. (35) and (37), we see that the point parti- 
cle source terms have been replaced by an integral over 
a worldtube surrounding the small body, as we desired. 
And the volume integral over the interior of does not 
diverge in f2, as it would in Eq. (37), because the region 
of integration excludes the interior of the worldtube. 

Finally, one should note the essential character of these 
integrals. They provide a type of Kirchoff representation 
[29, 65, 75] of a solution to the wave equation (135). 
However, while the integral representation is satisfied by 
any solution to the associated wave equation, it does not 
provide a solution. That is, one cannot prescribe arbi- 
trary boundary values on F and then arrive at a solution. 
The reason is that the worldtube is a timelike boundary, 
which means that field data on it can propagate forward 
in time and interfere with the data at a later time. How- 
ever, by applying the wave operator Eafi onto equation 
(140), we see that the Kirchoff representation of h is guar- 
anteed to satisfy the wave equation at each point x G il. 
In other words, the problem arises not in satisfying the 
wave equation in a pointwise sense, but in simultaneously 
satisfying the boundary conditions. However, since the 
tube is chosen to lie in the buffer region, these boundary 
conditions can be supplied by the buffer-region expan- 
sion. This can presumably be accomplished in a variety 
of ways, two of which I will discuss presently. Note that 
since the buffer-region expansion has been made to sat- 
isfy the Lorenz gauge to some order in using it as 
boundary data will enforce the Lorenz gauge in 12 to the 
same order. 

Now, recall that in almost all the derivations of the 
gravitational self- force (excluding those in Refs. [20, 24]), 
the first-order external perturbation was assumed to be 
that of a point particle. This was justified to some ex- 
tent by an argument first made by D'Eath [6, 10] and 
later used by Rosenthal [63]. The argument is based on 
the integral Eq. (138) and the asymptotically small size 
of the worldtube. First, note that the directed area ele- 
ment on the worldtube behaves as ^ ^"^{—n^ ). Also, in 
constructing the external solution, we formally assume 
r 1 (since the limit is constructed with fixed coordi- 
nate values in mind), which means that we can treat the 
Green's functions and its derivatives as quantities of or- 
der unity. Thus, the dominant term in the worldtube 
integral is determined by the derivative of the m/r term 

in ; using the result from the buffer-region expansion, 
this yields 



Hence, the boundary integral can be written as 



2m 

— {2Ua'Uf3' + ga'p') 



''Ea/3 



, (1) 



^ J mGap^ ^ {2Ua'Uf3' + Qa' f}')dt' d^' 



, (1) 



(142) 



where h-^^p is the contribution from the initial data sur- 
face E. Expanding the Green's function on the worldtube 
about the worldline 7, this becomes 



2mG^f,^p{2u'^u^ + g'''^)di 



(143) 



(141) 



where the barred coordinates correspond to points on the 
worldline, and i is proper time, running from t = to i ~ 
1/e. Equation (143) is the solution to the wave equation 
with a point particle source — except for the corrections 
of order It can be put in the more usual form of 
Eq. (35) by using the identity (A8). 

In the original derivation presented by D'Eath [6, 10], 
M was set to zero with no explicit justification. In Rosen- 
thal's later derivations [63], this step was justified based 
on the notion that we are interested in the limit in which 
the small body shrinks to a point. However, if the size 
of the body vanishes, then so too does its mass, in which 
case there is no perturbation at all; and at second or- 
der, setting 3£ to zero would create a divergent solution. 
Hence, discarding the order-^ corrections based on this 
argument is not justified. We could also argue that the 
order-^ terms must be discarded because the external 
solution cannot depend on the arbitrary radius of the 
tube. However, this second argument is also specious: 
One could just as easily express Eq. (142) as an inte- 
gral over any curve in the interior of F, rather than the 
central curve 7. But if one did so, then one would intro- 
duce mass dipole terms into the metric, and an explicit 
calculation of the error terms would show that they do 
not vanish. In some sense, this correctly implies that 
the choice of worldline at leading order is inconsequen- 
tial, since any choice within the worldtube results only 
in the introduction of a mass dipole, which is a second- 
order term, and the self-force will by definition set the 
resulting mass dipole to zero. However, this resolution 
becomes murky when we consider that the size of the 
tube must be left arbitrary to achieve a valid solution, 
and the mass dipole in the buffer region calculation is 
precisely order e^, rather than order e^. 

Instead, I present here an alternative argument to jus- 
tify D'Eath's conclusion: Suppose we take our buffer re- 
gion expansion of to be valid everywhere in the inte- 
rior of F (in M.E)^ rather than just in the buffer region. 
This is a meaningful supposition in a distributional sense, 
since the 1 /r singularity in is locally integrable even 
at 7. Note that the extension of the buffer-region ex- 
pansion is not intended to provide an accurate or mean- 
ingful approximation in the interior; it is used only as a 
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means of determining the field in the exterior. The rea- 
son I can do this is as follows: since the field values in 
CI are entirely determined by the field values on T, us- 
ing the buffer-region expansion in the entire interior of 
r leaves the field values in unaltered. Now, given the 
extension of the buffer-region expansion, it follows from 
Stokes' law that the integral over F in Eq. (138) can be re- 
placed by a volume integral over the interior of the tube, 
plus two surface integrals over the "caps" Jcap and Ecap, 
which fill the "holes" in J' and E, respectively, where 
they intersect F. Schematically, we can write Stokes' law 
as /i„t(r) = Ij^^, + - /r' ^here Int(F) is the inte- 
rior of F; this is valid as a distributional identity in this 
case.^^ The minus sign in front of the integral over F 
accounts for the fact that the directed surface element in 
Eq. (138) points into the tube. Because Jcap does not he 
in the past of any point in f2, it does not contribute to the 
perturbation at x € Hence, we can rewrite Eq. (138) 
as 

int(r) 
int(r) 

- C'/3'i^"'^'[Ga/3])dV^' + (144) 

where l^f!^^ is the contribution from the spatial surface 

t = EUScap. Now note that E"'^'[Gap] oc S{x, x'); since 
X ^ Int(F), this term integrates to zero. Next note that 

Ea'P'[fii^] vanishes everywhere except at 7. This means 
that the field at x can be written as 

C/3 = 1 G^/^'e^,0, [I^'JW + 4',^. (145) 

Int(r) 

Making use of the fact that Ea^lfi^] = 9^9c{l/r)}i^^p + 
0{r-'^), along with the identity d'^ddl/r) = -^■k5'^{x''), 
where 5'^ is a coordinate delta function in Fermi coordi- 
nates, wc arrive at the desired result 

C/3 = 2m / G^p^p{2u^u^ + g'''P)dt + h'^^^. (146) 
J J 

Thus, simply neglecting the O(^) terms in Eq. (143) 
yields the correct result, and in the region Q, the leading- 
order perturbation produced by the asymptotically small 
body is identical to the field produced by a point particle. 



^"^ Note that the "interior" here means the region bounded by F U 
Scup U Jcap- Int(r) does not refer to the set of interior points in 
the point-set defined by F. 



Gralla and Wald [20] have provided an alternative 
derivation of the same result, using distributional meth- 
ods to prove that the distributional source for the lin- 
earized Einstein equation must be that of a point parti- 
cle in order for the solution to diverge as 1/r. One can 
understand this by considering that the most divergent 
term in the linearized Einstein tensor is a Laplacian act- 
ing on the perturbation, and the Laplacian of 1/r is a 
flat-space delta function; the less divergent corrections 
are due to the curvature of the background, which dis- 
torts the flat-space distribution into a covariant curved- 
spacetime distribution. 

At second-order, the above method can be used to sim- 
plify Eq. (139) by replacing at least part of the integral 
over F with an integral over 7. I will not pursue this 
simplification here, however. Instead, I will present an 
alternative means of determining the metric perturba- 
tion. This method is based on a direct calculation of the 
boundary integral in Eq. (140). As such, it is somewhat 
similar in spirit to the Direct Integration of the Relaxed 
Field Equations (DIRE) used by Will et al. in post- 
Newtonian theory [76]. While the method used above 
relied on E^^y^^ being well defined as a distribution, a 
direct integration of the boundary integral can be per- 
formed, in principle, regardless of the behavior of h in 
the buffer region. Hence, it might be used at any order 
in perturbation theory. 

The method of direct integration proceeds as follows. 
As noted above, the Kirchoff representation of the so- 
lution is guaranteed to satisfy the wave equation at all 
points in 17, but it provides a valid solution only if, in 
addition, it agrees with the data on the boundary 9f2. 
Thus, the Kirchoff representation is guaranteed to be a 
solution in ^ if it satisfies the consistency conditions 

lim haji = ha' 13' 

for.T'eF. (147) 

lim n^\Iyiha(3 = \'ha'p' 

X — >x' 

However, these conditions allow h to contain a term such 
as (r — ^)^ ln(r — ^); both the term itself and its first 
derivative vanish in the limit r — > but the second- 
derivative does not. Since we seek a solution that is 
smooth and independent of I demand that h sat- 
isfy the following, stronger condition: Since the radius 
^% of the tube is small, the boundary data h' can be ex- 
pressed as an expansion in powers of ^ and £ — this is the 
buffer-region expansion. If x is near the worldtube, then 
r ^ meaning that hap can similarly be expressed as 
an expansion in powers of r and e. Defining <I>s(/(s)) to 
be an expansion of / for small s, I write the expansion 
of the boundary values as <i>e($^(/i')), and I write the 
expansion of the integral representation of the solution 
in as $e($r(/i))- I demand that these expansions are 
identical: 

$,($^(/i'))U^, = W)- (148) 

Hence, by expanding the integral representation of the 
perturbation near the worldtube and insisting that the 
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FIG. 4; The two-dimensional hypersurface S is defined by 
the intersection of the worldtube F with the past Ught cone 
of the point x. x is Unked to a point x' £ 5 by a null geodesic 
Of'. X and x' are separately linked to points x" = 7(f) and 
X = 7(t') by spacelike geodesies j3 and each of which is 
perpendicular to 7. 

result is consistent with the boundary data provided by 
the buffer-region expansion, all the free functions in the 
buffer region expansion will be determined. 

Since the equation of motion depends only on first- 
order terms, for the purposes of this paper I will limit the 
expansion just described to first order. The expansion is 
performed only in the buffer region, meaning that it pro- 
vides an explicit expression for the perturbation only in 
that region. However, by imposing the consistency con- 
dition, the boundary data on T can be determined to any 
desired order of accuracy in M] using this boundary data, 
the solution in 57 will then be determined to the same or- 
der of accuracy. A similar procedure could be adopted at 
second order and above. At those orders, the expansion 
of the boundary integral would yield ^-dependent terms 
that would be grouped with the volume integral over 17; 
this combination would yield an approximation to the in- 
homogenous part of the solution. The homogenous part 
of the solution would be dealt with in the same manner 
as the first-order perturbation. 



B. The boundary integral 

Since the calculation in this section is intended primar- 
ily as a proof of principle, rather than calculating lij^ 



I will calculate its approximation /i' ' ; in other words, 
I will consistently neglect acceleration terms. Hence I 
take the boundary data on the tube to be defined by 
h = -ii^' ^' -I- eft,*^'"' + erli'"''\ and the field outside the 
tube to be the expansion of 

an 

+ 0{e^) (149) 

to order er. Since the volume integral contributes only 
0[e^,r^) terms, it is neglected here. 

Two parts of the boundary lie within the causal past 
of x: the spatial hypersurface E and the worldtube F. 
The contribution to the field from the data on E is given 
by 

^Sal3 = ^ j [Gaff ^ flj'S'-f,' " ^V'^'^a/ ^ ' 

(150) 

where the data h^,g, is constrained to satisfy the Lorenz 
gauge and merge smoothly with the buffer region expan- 
sion. I assume that Ii^! can be expanded in a regular 
power series in r, 

m>0 

and that each ht^^j^ can be decomposed into irreducible 
STF pieces. Because this data can only contribute to the 
homogenous, free functions in the buffer region expan- 
sion, we can infer the nonzero pieces of h-^^^ from that 
expansion: 

(152) 

C=^r, (153) 
C: = <5,,jr^^'°'+^4i->, (154) 
4^=<-^>n\ (155) 
4a = + '^^^r/ + ea^&^fn^, (156) 

Now consider the integration over F. In the buffer 
region, I define A(e) such that r ~ ^ ~ ^(^)- Since 
the function A(e) is arbitrary, except that it must van- 
ish in the limit e ^ 0, we can use e and A as in- 
dependent expansion parameters. The volume element 
on F is given by dS*^' = -nf,'N{x')^'^dt'dn' , where 
N{x) = 1 + ^£cdit)x'"^ + 0(A3,e), and t', ^, and 6*'^ 
are Fermi coordinates based at 7. The boundary data is 
constructed from 

(1) 1 (1,0) _ (1.1) / .-.x 

S'V = ^^'5' + + ^'^-r'S' + 0{X') (158) 
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where /ly^?, fi^'°s', and /ly^, are obtained by setting the 
acceleration to zero in Eqs. (80), (82), and (91). 

The integral over F can be divided into two regions: the 
convex normal neighbourhood Af of x — consisting of all 
the points that are connected to x by unique geodesies — 
and the complement of the convex normal neighbour- 
hood. In Af, the Green's function admits the Hadamard 
decomposition [29] 

G„/'' = U^/''S+{a) + Vo,i3'''''0+{-a), (159) 

where a{x,x') is Synge's world function, which is equal 
to one-half the squared geodesic distance between x and 
x' . Derivatives of this biscalar will be denoted by, e.g., 
CT^ = cr.^. The delta function S^{a{x,x')) has support 
on the past light cone of x, while the Heavisidc function 
9+{—a{x,x')) has support within the past light cone. 

Substituting these expressions into the boundary inte- 
gral, we find that it can be broken into several pieces; 
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+ J t)ijjNdt'dn' + + 0{e') , (160) 

(r\A^)n7" 

where / is the past of x, and 6' is the derivative of the 
delta function. Inside the normal neighbourhood, the 
terms in the integrand are given by 
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+ 0(A5), (163) 

where V„' = n'^Va'. The "direct" and "tail" titles should 
be self-explanatory. Outside the normal neighbourhood, 
the term in the integrand is 



in /i* ' ; each of the above expansions is performed to 
an order sufhcient to meet this requirement, given that 
(5+(ct) - 1/A2 and 8\{a) - l/A^. 

It is convenient to adopt cr as an integration variable, 

which can be done using the transformation dt' = — , 

where CTa' (x, x')it" can be thought of as a measure 
of the luminosity distance from x to x' . Note that the 
four-velocity at x' is defined as the tangent to a curve of 

constant r and Q^: that is, u" = ^7-|p. Since t' is the 
proper time on the worldline, rather than the proper time 
on the generators of the worldtube, this four-velocity is 
not normalized. (This implies that ^ is not an affine 
parameter on the geodesic connecting x to a;'.) 

After performing this change of variables, we eliminate 
the (5' term in the boundary integral by using the identity 



\)tfT^i\<j)di!d^ 



Here S is the intersection of the past light cone of x 
with the worldtube. For simplicity, I assume that the 
normal neighbourhood of x is large enough for S to be 
well defined and for the intersection 5 n A/" to be empty. 
I also assume that x is late enough in time for S to be 
closed, such that it has the topology of a sphere. (If x 
is not sufficiently late in time, then S will be "cut off' 
where it intersects E.) 
We can now express h as 
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where 
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The metric perturbation has three types of contributions: 
the "direct" type arising from data on the light cone; the 
"tail" part arising from the interior of the light cone; and 
the contribution from the initial data surface. 

The "direct" and "tail" contributions will be calculated 
explicitly in the following subsections. We begin with the 
direct contribution. 



r-tail 
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These expressions are completely general; they can be 
simplified by making use of the fact that Vi'/iy^V = 
0{\\Xe) = V„./i';;";. Note that , \)%\ and \)tf 
are scalars at x' and rank- two tensors at x. Also note 
that we require the final answer to be accurate up to er- 
rors of O(A^), in order to determine the free functions 



C. Integral over the past light cone 

Each of the quantities in the bitensor f)^'^ can be 
expanded in powers of A by first expanding the x' - 
dependence about the point x = 7(i') and then expand- 
ing the x-dependence about the point x" = 7(t). See 
Fig. 4 for a depiction of the relationship between these 
points. The expansions are provided in Appendix C. 
Most significantly, the distance and its time-derivative 
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dt' i- are expanded as ^ = 
and ^ = i^^Xix^y? . 
terms are the flat-spacetime values 



A(^?o + A^i + A2^ + ...) 

, where the leading-order 



(168) 
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After making use of these two expansions, and express- 
ing }i (t') and }i ' ' (t') in terms of their values at we 
can express i)^" explicitly as 
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^0 ^ ^0 



Each term in this expression is further expanded using the results of Appendix C, which details the expansion of the 
Green's function and cr^'(a;, x'). In order to integrate the final, fully expanded expression for t)*^" over the surface S, 
we make use of the angular integrals displayed in Appendix B 1 . 



The end results of these calculations arc as follows: The h^^ terms in [) contribute 
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The hi^p terms integrate to zero at the orders of interest. And the h^^J terms contribute 
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Note that the h^p terms are all ^-dependent, and they the tetrad (e°,e^) in terms of the coordinate one-forms 
are necessary to cancel ^-dependent terms arising from {ta,x^), we can write this result in Fermi coordinates as 
h^^p ' . All the actual terms that appear in the buffer 
region 

perturbatiuii, uuu tuc icguim tcimo, ouv^n ao 't^^ 

required on the boundary to ensure the consistency of 
the solution. 

Putting these results together, we arrive at 
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0(A2,£), 



(173) Note that the l/r term in this result agrees with the 1/r 
term that was used for boundary data. 



where all tensors are evaluated at time t. By expressing We now proceed to the calculation of the tail terms. 
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D. Integral over the interior of the past light cone I) can be written as 



The interior of the Ught cone covers the worldtube from 

the lower hmit t' ~ 0, where the worldtube intersects 

the surface E, to the upper limit defined by the surface 

S. Because the Fermi time varies over 5, I split the 

integration into two regions: one region from t' = Q to 

t' = t^in = min{t'}, and another from t' = t^^hi to t' = 
s 

tmax = max{t'}. The integral over the worldtube is then 
s 

expressed as 



dfl'dt' = dt' dn' 
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+ 0(a2) (179) 
0{\^). (180) 



In addition, in the second and third integral, which lie 
within the normal neighbourhood of x, the Green's func- 
tion can be replaced with V^^ " ^ and we can use the 
near-coincidence expansion 



^^ej^f\ej'RK\{t)+OiX) (181) 



In the integral from tmin to imax, the angles 9'^ = {9\ cj)') 
are cut off at some maximum values defined by S. 

Because f)*^'' is of order A", and we only seek terms 
up to order A, we can further simplify the integral. For 
x' G S, we can write the time difference t — t' as, t — 
t' = ^ + 0(A2), where the 0{\^) error term consists 
of acceleration and curvature terms (see Eq. (C18)). I 
choose 6' to be the angle between x"^ and a;", such that 
t' ^ t - -1-^2 _ 2r^ cos 61' + 0{\'^). From this we 
infer that the maximum and minimum times on S are 
given by tmax = < - (r - ^) + ©(A^) and tmin = t - {r + 
+ O(A^). The value for tmax corresponds to the time 
at the point on S closest to x; the value at tmiu is the 
time at the point furthest from x. Since the maximum 
value of 9' at a given value of t' is determined by the 
intersection with 5, it is given by 



cos 9n 
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Since <f>' runs from to 27r everywhere on S, its maximum 
value is 27r, independent of t' and 9' . 

Making use of these approximations, we can expand 
the first integral, running from t' = Q to t' = imin, about 
the upper limit t' = t. This enables us to write the 
integral over the worldtube as 



■t /'47r /'47r 

dn'dt' = j dt' I dn' -{r + .^) / dn' 

"'0 "'0 



+ d(t)' dt' / dcose'' + Ci(A^) . 

(177) 



rnr 



for x' near x" = 7(t). 

Substituting these expressions into the integral, and 
noting that / n'°-dVL' = 0, yields the result 

^ fCflNdn'dt'= I 2m(G„^oo + G„;3 di' 
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Note that the integrals are cut off at = < — 0+ to 
avoid the singular behavior of the Green's function at 
coincidence. 

Now, since x is near the point x" on the worldline, we 
can expand the integrand as 



= (Ga„p„g^p + Vy//G„„^„a^e2' : 
X gf gf (2u''u^ + g^f") (183) 

Substituting this into the integral results in the expansion 
^ / l)':^'Ndn'dt' ^ gfgf (/irt"0" + h^l.p.^.efx^ 
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where I have defined 



With the simplification Vn'hy^, 
X^i'/iy^/, along with the expansion 
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(186) 
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By making use of the identity (A8), we can express these 
tail terms in their more usual form: 



tail 



\gc."p''G'"s„^p)u^u''dl (187) 
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tail 



The complete tail term will consist of the sum of the /ip 
terms and the terms. 

In Fermi coordinates, the final result of this section is 

tail tail 
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+ 0(A2,e). (189) 

Note that the .^-dependent term in this equation exactly 
cancels the ^-dependent term in Eq. (174). In addi- 
tion, note that this expansion is identical to the one in 
Sec. IV A only after explicit factors of the acceleration 
are set to zero. This means, in effect, that when compar- 
ing individual components of our expansion here to those 
in our previous expansion in the buffer region, we should 
replace the covariant derivate in the Fermi-coordinate ex- 

tail 

pression for h^-^n ^n^n with a partial derivative. 



E. Identification of unknown functions 

We now combine the results of Eqs. (150), (174), and 
(189) to arrive at an expansion of the form 
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which we will identify with the expansion defined by 
Eqs. (80), (82), and (91). After defining the tail terms 



tail 


tail 






tail 


tail 




— '^TIJc 



, (1,0) 



(1,1) 



(191) 
(192) 



and decomposing the results into STF pieces, we find 
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TABLE II: Symmetric trace-free tensors in the first-order 
metric perturbation in the buffer region, written in terms of 
the electric-type tidal field £ab and the tail of the perturba- 
tion. 
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mSab 
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After setting explicit factors of the acceleration to zero, 
this expansion agrees with Eqs. (80), (82). and (91). By 
comparing the two sets of equations, we identify all the 
unknown STF tensors in the buffer region expansion. 
The results of this identification are listed in Table II. 
Note that these identifications are modulo the accelera- 
tion that appears in the covariant derivative in h^^p^- 

Note that for this solution to agree with the results of 
the buffer region expansion, it must satisfy the relation- 
ships given in Eqs. (89) and (90). In terms of the tail 
integral, these relationships read 
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tail 



(198) 
(199) 



where it is understood that the equations hold only for 
a = 0. By using the Green's functions identities (A6), 
(A7), and (A8), and neglecting acceleration terms, one 

can easily show that the tail terms hj, satisfy these re- 
lationships. Hence, we must constrain the initial data 
terms lij! to independently satisfy them, which implies 

^dt + sje-^'-"') , (200) 

(201) 
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The reader should take note of two important facts 
about the metric perturbation derived here. First, as we 
expected, the expansion displayed above is identical to 
the expansion of the point particle solution in the neigh- 
bourhood of the worldline. Second, and again as we ex- 
pected, the expansion is completely determined by the 
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most singular, e/r, term in the metric. Although the 
nonsingular terms are required to maintain consistency 
at the boundary, one can derive all of them simply by 
using the 1/r term as boundary data. 

Now, the principal purpose of the calculation of the 
boundary integral was to express the equations of motion 
in terms of the body's past history. The correction to the 
body's mass, given in Eq. (120), can now be written as 

tail tail 

Sm{t) = Sm{0) + im/iQo + f^mS''''h^^ . (202) 
In covariant form, this is 

Sm{t) = Sm{0) + (Sg"'' + h^p. (203) 

This is similar, but not identical to, a result found by 
Mino, Sasaki, and Tanaka [18]. The source of disagree- 
ment between the two results is not clear. It is worth 
noting that both results appear at one order lower than 
that given by Thorne and Hartle [59] , who chose to elim- 
inate the homogeneous field . 

The leading-order acceleration of the body, given in 
Eq. (121), is 

a^^'^^Cl- Vo-;^^.Sa- (204) 

(Here, again, the right-hand side of this equation is to be 
evaluated at a = a*"' ~ 0.) In covariant form, this result 
can be written as 

a" = -hj {9'" + u-u') (2^ - h';^s) u^u^ 

+ :^R''0jSU^S-'' + 0{e^) (205) 
2m 

where I have again used the definition S'^^ = eje^^e'^'^^ Sj. 
The spin term is the usual Papapetrou spin force. The 
tail term is the usual MiSaTaQuWa self-force — except 
that the tail integral is defined as the sum of an integral 
over the worldline, cut off at i = 0, and an integral over 
an initial data surface. Of course, Eqs. (202) and (204) 
hold only in the Lorenz gauge. 

This concludes what might seem to be the most egre- 
giously lengthy derivation of the self-force yet performed. 
It is hoped, however, that along with the additional 
length has come additional insight. 

VI. DISCUSSION 

A. Summary 

In this paper I have presented a new derivation of the 
gravitational self-force for a small body. It is based on 
the familiar technique of using two expansions of the 
metric: an inner expansion that is more accurate near 
the body, and an outer expansion that is more accurate 
far from the body. However, unlike in earlier deriva- 
tions, I have formulated these expansions in terms of 



a fixed worldline 7 defined in the external background 
spacetime. The self-consistent equation of motion of this 
worldline then follows directly from solving the Einstein 
equation. When combined with the first-order metric 
perturbation, the equation of motion defines a solution 
to the Einstein equation accurate up to order errors 
over times t < 1/e. When combined with the second- 
order perturbation, it defines a solution accurate up to 
order errors on the shorter timescale ^ 1. 

My approach began with a general analysis of the Ein- 
stein equation, up to second order in the body's mass, in 
a buffer region around the small body. Since the buffer 
region is assumed to be free of matter, my calculation 
is valid only for bodies that are sufficiently compact to 
avoid tidal disruption. An equation of motion for the 
body's worldline was derived from the condition that the 
body must possess no mass dipole in coordinates cen- 
tered on the worldline. From this purely local-in-space 
analysis, we found an expression for the acceleration in 
terms of irreducible pieces of a homogeneous solution to 
the wave equation — the Detweiler- Whiting regular field, 
which is regular on the worldline. This homogeneous, 
regular field was not determined by the buffer-region ex- 
pansion, since it can be determined only by boundary 
conditions. 

A formal expression for the metric perturbation was 
obtained by casting the Einstein equation in a relaxed 
form, via the imposition of the Lorenz gauge. This re- 
laxed form can be solved iteratively, with the perturba- 
tion at each order given by the sum of (1) an integral 
over a region outside the body and (2) an integral over 
an initial data surface and a worldtube surrounding the 
body. Boundary data on the worldtube are provided 
by the buffer-region expansion. At first order, it can 
be shown that the integral representation is identically 
equal to the perturbation produced by a point particle 
moving on 7. At higher orders, because of the increasing 
singularity of the metric perturbation, only parts of it 
can be simplified in the same way. Because of this lim- 
itation, I introduced a method of direct integration. In 
this method, the Detweiler- Whiting regular field in the 
neighbourhood of the body is determined, in terms of ini- 
tial conditions and the body's past history, by expanding 
the integral representation in the buffer region and de- 
manding its consistency with the boundary data on the 
tube. 

An essential assumption in this derivation is that the 
acceleration of the fixed, e-dependent worldline possesses 
an asymptotic expansion beginning in powers of s. In ad- 
dition, I made the following assumptions: the exact met- 
ric possesses asymptotic expansions of the form given in 
Sec. IIIC, there is a smooth coordinate transformation 
between some internal local coordinates and the exter- 
nal Fermi coordinates in a neighbourhood of the world- 
tube, the Lorenz gauge condition can be imposed ev- 
erywhere in the region of interest, and the expansion of 
the metric perturbation satisfies both the wave equations 
and (when combined with the expansion of the acceler- 
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ation) the gauge condition at fixed functional values of 
the worldline z^{t). While these, especially the last, are 
strong assumptions, they undoubtedly lead to an emi- 
nently useful, systematic approximation scheme. It is 
worth repeating that while the choice of gauge is not es- 
sential in finding an expression for the force in terms of 
the field in the buffer-region expansion, it is essential in 
my method of determining the field itself. Without mak- 
ing use of the relaxed Einstein equations, no clear method 
of globally solving the Einstein equation presents itself. 

One fruitful avenue of further research might be to 
explore methods of solving the Einstein equation in al- 
ternative gauges but still within the context of a singular 
expansion that holds 7 fixed. This might require further 
thought on the behavior of gauge transformations in such 
an expansion. However, such details of the formalism are 
most likely to be made sense of not at the level of the 
field equations, but at the level of the action. Since par- 
tial derivatives do not act directly on the worldline, the 
dependence of the metric perturbations on the worldline 
does not appear directly in the field equations. But at 
the level of the action, using functional derivatives, the 
role of the worldline becomes transparent. This will be 
explored in a future publication [43]. 

It is also worth noting that the methods used here 
would work in many other cases. For example, the direct 
calculation of the boundary integral can be used to com- 
pletely determine the force even if the source cannot be 
represented as a distribution. Also, these methods could 
be used to derive self-consistent equations of motion for a 
charged body; the expansion of the acceleration in powers 
of £ would automatically yield an order-reduced equation 
of motion, with no runaway solutions (c.f. the recent cal- 
culation by Gralla, Harte, and Wald [77]). 



B. Comparison with alternative methods 

One of the goals of this paper was to construct 
an approximation scheme that closely mirrors the ex- 
tremely successful methods of post-Newtonian theory 
[3, 4, 55, 56, 62, 64, 76]. As such, many of the methods 
used here are similar to those used in post-Newtonian ex- 
pansions. For example, the expansion with a fixed world- 
line meshes well with the use of the relaxed form of the 
Einstein equation [45, 46], which can be solved without 
specifying the motion of the source, and which is the 
starting point for post-Minkowski and post-Newtonian 
expansions. And the use of an inner limit near the body 
corresponds to the use of the "strong-field point parti- 
cle limit" used by Futamase [55, 56]. In addition, the 
calculation of the motion of the body in this paper is 
somewhat similar to the methods used by Futamase and 
others [3, 55, 56, 62, 64], in that it is based on a multipole- 
expansion of the body's metric in the buffer region. Fi- 
nally, the direct integration of the relaxed Einstein equa- 
tion mirrors the approach of Will et al. [76] . 

There are, of course, differences between the two cases. 



In particular, when the finite size of the body is taken 
into account in post-Newtonian theory, because the back- 
ground is fiat, finding an equation of motion for the mass 
dipole of the body is equivalent to finding an equation of 
motion for its worldline. Although this method, or meth- 
ods similar to it, has also been used in curved spacetimes 
[24, 25, 59], it is somewhat problematic because the mass 
dipole corresponds to a displacement from the center of 
a given coordinate system. But in a curved spacetime, 
such a displacement is meaningful only when it is in- 
finitesimal. Of course, if at a given instant the coordinate 
system is mass-centered, then the second time-derivative 
of the mass dipole is equivalent to the acceleration of 
the worldline; but since there is no unique global time 
in a curved spacetime, it is more meaningful to speak of 
a curve about which the body is centered for its entire 
history, rather than just at a given time. 

The more significant goal of this paper was to develop 
a unified and self-consistent formalism to treat the grav- 
itational self-force problem. Because the problem con- 
sists of solving singular perturbation equations, I have 
emphasized the foundation of the formalism in singular 
perturbation theory. Because the formalism uses a self- 
consistent worldline and a finite sized body, it is (poten- 
tially) valid on both short and long timescales, and both 
very near to and far from the small body. As such, it can 
be used to study (or incorporate studies of) the space- 
time near the small body, the long-term motion of the 
body, and the perturbations produced by it, including 
the gravitational waves emitted to infinity. 

This contrasts with the most recent derivation of the 
self- force, performed by Gralla and Wald [20] . In terms of 
the concrete calculation of the force in the buffer region, 
my calculation is very similar to theirs, though it differs 
in many details. (One such difference is that the pertur- 
bation I derive satisfies the Lorenz gauge at all orders 
in r in the local expansion, whereas Gralla and Wald 
do not impose the Lorenz gauge on the most singular, 
order-£^/r^, term in their calculation.) However, their 
approach constructs a regular expansion in which both 
the worldline and the metric perturbation are expanded; 
they suggest that in order to arrive at a self-consistent 
set of equations, one must make a "leap of faith" from 
the results of their regular expansion. I instead take the 
stance that the self-consistent equation of motion can, 
and should, be justified by a more systematic approach; 
and I have presented one such approach in this paper. 
From the results of this approach, one can easily derive 
the results of the regular expansion: simply by expand- 
ing the £-dependent worldline, one derives a leading-order 
metric perturbation sourced by a particle on a geodesic 
(plus secularly growing corrections); and the usual steps 
involved in deriving the geodesic-deviation equation leads 
to an equation of motion for the deviation vector "con- 
necting" the geodesic to the exact worldline. Contrari- 
wise, one cannot derive the results of the singular expan- 
sion from those of the regular expansion. 

Other methods have been developed (or suggested) 
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to accomplish the same goals as my own. One such 
method is the two-timescale expansion suggested by Hin- 
derer and Flanagan [35]. As discussed in Sec. II C, their 
method continuously transitions between regular expan- 
sion, resulting in a global, uniform-in-time approxima- 
tion.^^ This method contrasts with the one presented 
in this paper, in which a single expansion has been con- 
structed by treating the worldline as fixed. 

The fundamental difference between the two methods 
is the following: In the two-timescale method, the Ein- 
stein equation, coupled to the equation of motion of the 
small body, is reduced to a dynamical system that can 
be evolved in time. The true worldline of the body then 
emerges from the evolution of this system. In the method 
presented here, I have instead sought global, formal solu- 
tions to the Einstein equation, written in terms of global 
integrals; to accomplish this, I have treated the worldline 
of the body as a fixed structure in the external spacetime. 
The two timescale method is, perhaps, more practical 
for concrete calculations, while the global solutions pre- 
sented here are primarily of formal interest. However, 
the two methods should agree. Note, though, that Hin- 
dercr and Flanagan have identified transient resonances 
in EMRI systems, which lead to half-integer powers of 
£ in their asymptotic expansions. It is not clear that 
such effects are correctly accounted for in the method 
presented in this paper. 



C. Prospects for a global solution 

The principal practical goal of solving the self-force 
problem is to find the waveform emitted from an EMRI. 
In order to extract the parameters of an EMRI system 
from its waveform, we must have a model that tracks the 
wave's phase to within an error of order e over a time 
period 1/e. This presents several problems. 

First amongst these problems is the potential for sec- 
ular errors. For example, secular errors might arise due 
to ignoring the slow evolution of the background space- 
time. Throughout this paper, I have assumed that the 
external background metric is e-independent. However, 
in practice, it might possess a slow time dependence that 
would account for the backreaction of the perturbations 
on the background spacetime; for example, in an EMRI, 
the large black hole's absorption of gravitational waves 
slowly alters its mass and spin parameters. Any such 
effect leaves the expression for the self-force unchanged, 
and it can be easily incorporated into the perturbations 



Note that simply patching together a sequence of regular ex- 
pansions, by shifting to a new geodesic every so often using the 
deviation vector, would not accomplish this: such a procedure 
would accumulate a secular error proportional to the number of 
"shifts" multiplied by a nonlinear factor depending on the time 
between "shifts." And this error would, formally at least, be of 
the same magnitude as the solution itself. 



presented here [43]. However, an equation for the slow 
evolution itself is unknown. Presumably, it can be de- 
termined from an averaged version of the Einstein equa- 
tion, of the form {E^^[h]) = 2{R^^) +2{5'^R^^) + .... 
In an EMRI system, the average of the wave operator 
will most likely vanish, because the body's orbit is quasi- 
periodic. The averaged equation will then relate 5^R 
to the background Ricci tensor i?, as in the pioneering 
work of Isaacson [78]; this corresponds to the effect of 
quadrupole radiation on the background. In practice, 
the averaged equation might be solved by using some 
ansatz for the background metric — e.g., the Kerr met- 
ric with slowly varying mass and spin parameters. The 
feasibility of such a calculation is unclear; the need to 
perform it will most likely be determined by examining 
the magnitude of secular growth in a solution that ignores 
backreaction. See Refs. [35, 79] for more information on 
the backreaction in the self-force problem. 

Putting aside the backreaction problem, other secular 
errors will also arise due to neglected terms in the accel- 
eration and metric perturbation. Although the approach 
taken in this paper is designed to avoid such errors, a 
concrete implementation will nevertheless contain them. 
I have defined the worldline as a fixed curve; proceed- 
ing to successively higher orders in perturbation theory 
yields successively more accurate equations of motion for 
this curve. However, if we stop at any given order and use 
any given equation of motion, then the worldline based 
on that equation of motion will deviate secularly from 
the true worldline. This in turn implies that the metric 
perturbation will accumulate secular errors. 

Hence, we must have an equation of motion that lim- 
its these errors to 0(e) after a time \/e. If we use the 
first-order equation of motion, we will be neglecting an 
acceleration ~ e^, which will lead to secular errors of 
order unity after a time Thus, the second-order self- 
force is required in order to obtain a sufficiently accurate 
waveform template. In order to achieve the correct 
waveform, we must also obtain the second-order part of 
the metric perturbation; this can be easily done, at least 
formally, using the global integral representations outside 
a worldtube. 

A formal expression for the second-order force has 
already been derived by Rosenthal [21, 63]. However, 
he expresses the second-order force in a very particular 
gauge in which the first-order self-force vanishes. This is 
sensible on short time scales, but not on long timescales. 
Furthermore, it is not a convenient gauge, since it does 
not provide what we wish it to: a correction to the 
nonzero leading-order force in the Lorenz gauge. 

Thus, we wish to obtain an alternative to Rosenthal's 
derivation. Based on the methods developed in this pa- 



Proceeding to second order will also be useful for examining other 
systems, such as intermediate mass ratio binaries, over shorter 
timescales. 
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per, there is a clear route to deriving the second-order 
force. We would construct a buffer-region expansion ac- 
curate up to order e^. Since we would require the order 
e^r terms in this expansion, in order to determine the 
acceleration, we would need to increase the order of the 
expansion in r as well: specifically, we would need terms 
up to orders e^r^, er^, e'^r, and e^r^. These could be 
calculated using the methods presented previously. 

Unfortunately, such a procedure could be prohibitively 
difficult. Hence, we might consider a much simpler alter- 
native: the method of matched asymptotic expansions. 
Using this method, we would need the buffer region ex- 
pansion to be accurate to order £^r — the order at which 
the second-order acceleration appears in the background 
metric — meaning that we would need to extend the buffer 
region expansion by one order in r, but not in e. An in- 
ternal solution is already known to that order [58] , mean- 
ing that only the external solution must be found. The 
equation of motion would then be determined by find- 
ing a unique coordinate transformation that makes the 
external and internal solutions identical in the buffer re- 
gion. As will be discussed in Ref. [43], this method is 
somewhat problematic, because the coordinate transfor- 
mation is not unique. However, it should be possible to 
overcome this problem, and a calculation of the second- 
order force by this means is entirely feasible. 

Even if we can obtain an aproximation with the desired 
accuracy on the timescale 1/e, there remains at least one 
additional difficulty. The waveform itself is to be calcu- 
lated at future null infinity, J^^. At first glance, it might 
seem that we can extend the size of our domain such 
that its future null boundary J' is pushed out to J^~^ at 
one end and to the event horizon of the large black hole 
at the other. However, the size of our domain is intended 
to be of size 1/e. Thus, if we wish to enlarge the domain 
out to infinity, we must match the solution within it to 
an outgoing wave solution at its future null boundaries. 



D. Conclusion 

Throughout this paper, I have taken the stance that 
finding a useful approximate solution to the exact Ein- 
stein equation, such as that provided by singular per- 
turbation theory, is more important than finding an ex- 
act solution to the approximate Einstein equation, such 
as that provided by regular perturbation theory. In the 
gravitational self-force problem, a useful approximate so- 
lution is one that remains valid on long timescales, self- 
consistently incorporates the acceleration of the small 
body, and accounts for its asymptotically small, but fi- 
nite, size. In this paper, I have developed a formal ap- 
proximation scheme that promises to satisfy these crite- 
ria, and which can be extended to any order in pertur- 
bation theory. 

However, I have also taken the stance that an approx- 
imate solution of the exact Einstein equation must be 
an approximation to an exact solution if it is to ren- 



der a meaningful test of General Relativity. As such, I 
have emphasized how the singular expansions developed 
in this paper might be related to an exact solution. A 
far more rigorous, technical, and perhaps altogether un- 
feasible study would be required to show whether or not 
the asymptotic solution developed here actually does ap- 
proximate an exact solution. 

Of course, even if the solutions are proven to be asymp- 
totic approximations, they remain purely formal. A prac- 
tical calculation of the motion of a small body will most 
likely require a numerical implementation, which will re- 
quire a formulation of the wave equation, coupled to an 
equation of motion for the source, that is viable for nu- 
merical calculations. 
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APPENDIX A: GREEN'S FUNCTIONS 

I follow the notation and conventions of Ref. [29] . The 
Green's function for the tensor wave operator is defined 
by the equation 

(Al) 

where S{x, x') = 5'^{x°'—x'^ )/ \/\g\, and is the parallel 
propagator from x to x' . The Green's function for the 
vector wave operator is defined by 

{g'l.U - i?/) G/ {x, x') = -ing^:;S{x, x'), (A2) 

and the Green's function for the scalar wave operator is 
defined by 

{D^ XR)G{x,x') ^ -AnSix^x'). (A3) 

All quantities are defined with respect to the background 
metric g. 

If the point x' lies in the convex normal neighbourhood 
of X, then the retarded gravitational Green's functions 
can be written in the Hadamard decomposition 

GaPa'p' — Uaf3a'l3'6+{(T) + Vafja' 13' +{ — (t) , (A4) 

where (ct) is a delta function with support on the past 
light cone of x, 0^{—a) is a Hcavisidc function with sup- 
port in the interior of the past light cone of x, a is Synge's 
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world function, equal to one-half the squared geodesic 
distance between x and x', and Uapa'p' and Vapa'p' are 
smooth bitensors. The advanced Green's function has 
an analogous decomposition, with 5+{a) and 9+{—a) re- 
placed by 5 -{a) and 9 -{—a), which, respectively, have 
support on and within the future lightcone of x. The sin- 
gular Green's function, which satisfies the same defining 
equation as the retarded and advanced Green's functions 
but which has support only on and outside the past and 
future light cones, is given by 



If g is a vacuum metric, such that R = = R^^^ 
one can easily derive the following identities: 



(A5) 
then 

(A6) 
(A7) 
(A8) 



Equation (A6) follows from taking the divergence of 
Eq. (Al) and the covariant derivative of Eq. (A2). Equa- 
tion (A7) follows from taking the divergence of Eq. (A2) 
and the covariant derivative of Eq. (A3). Equation (A8) 
follows from contracting the primed indices in Eq. (Al). 
In each case, these operations show that the two relevant 
bitensors satisfy the same differential equation; the equa- 
tions (A6)-(A8) hold when the bitensors on the left and 
right satisfy identical boundary conditions. 

Equation (A7) appears in Ref. [37]. To the best of my 
knowledge, Eqs. (A6) and (AS) have not been presented 
in the literature. 



be found from the formula 
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These coefficients can then be decomposed into irre- 
ducible STF tensors. For example, for s = 1, we have 



(+) 



(-) 



where the T*^"^ 's are STF tensors given by 
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Similarly, for a symmetric tensor Ts with s = 2. we have 



T, 



ab(L) 



T. 
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where 
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APPENDIX B: STF TENSORS 

This appendix briefly reviews the use of STF decom- 
positions and collects several useful formulas. Refer to 
Ref. [66, 80, 81] for thorough reviews. All formulas in 
this section are either taken directly from Refs. [80] and 
[66] or are easily derivable from formulas therein. 

Any Cartesian tensor field depending on two angles 6^ 
spanning a sphere can be expanded in a unique decom- 
position in terms of symmetric trace-free tensors. Such a 
decomposition is equivalent to a decomposition in terms 
of tensorial harmonics, but it is sometimes more conve- 
nient. It begins with the fact that the angular depen- 
dence of a Cartesian tensor Ts{9^) can be expanded in 
a series of the form 



(Bl) 



e>o 



where S and L denote multi-indices S = ii-.-ig and 
L = angular brackets denote an STF combination 

of indices, n° is a Cartesian unit vector, = n^^ . . . n^* , 
and = n^^\ This is entirely equivalent to an expan- 
sion in spherical harmonics. Each coefficient Tgi^^ can 



T 



(+2) 



T 



L+2 

( + 1) = 
L+l 



T, 



T 



L-l 

(-2) 
L-2 

Kl 



(L+2), 
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STF(r(p,^,)gi_iei,^/«), 



(0) _ 



£ + 2 L+i 

U{21 ~ 1 
(£ + l)(2^-f 3) "I 



(B8) 
(B9) 

(BIO) 



_ 2{l~\){2l^\) 



(£+l)(2^+l) STF(ro-rt,^L-2e.._/'^) 
2£- 3 



2(. + l 
3 3'^- 



(Bll) 
(B12) 
(B13) 



These decompositions are equivalent to the formulas for 
addition of angular momenta, J = S + L, which results 
in terms with angular momentum i — s < j < £ + s; the 
superscript labels (±n) in these formulas indicate by how 
much each term's angular momentum differs from £. 

By substituting Eqs. (B3) and (B7) into Eq. (Bl), we 
find that a scalar, a Cartesian 3- vector, and the symmet- 
ric part of a rank-2 Cartesian 3-tensor can be decomposed 
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as, respectively, 



(B14) 



i>0 



Ta{0^)=Y,BLnaL 



L-l 



71^ ^ +e'aiD,L-in- 



1>1 



(B15) 



e>o 



e>i 



{aJ-b)jL-2ni 



(B16) 



Each term in these decompositions is algebraicaUy inde- 
pendent of aU the other terms. 

We can also reverse a decomposition to "peel" a fixed 
index from an STF expression: 

+ 1) STF T,(L) = T^L) + e STF T,^i^,l-i) 



2£ 



21- 



-STYT\JL-l)K^■ (B17) 

L 1^ 



In evaluating the action of the wave operator on a de- 
composed tensor, the following formulas are useful: 



n n = n 



2e + i 



Ucrf^ = — n^, 

2e + i 



(B18) 
(B19) 



lU + 1) 

rdcfiL = ~£hcL + . , 5c{hni2...ii) , (B20) 



21 ^\ 



-n , 



'^'^^^ ~ (2£+l) " ■ 



(B21) 
(B22) 
(B23) 



In evaluating the i-component of the Lorenz gauge con- 
dition, the following formula is useful for finding the most 
divergent term (in an expansion in r): 



^ 2£ + 1 ^ 

£>2 



(B24) 



And in evaluating the a-component, the following for- 
mula is useful for the same purpose: 



ah 



(^ + 2)(^ + 3) , ' 

2i + 3 ^ 6 L 



E 



.>-t L2(2^ 

+E 



-L-l 



£>1 



(^ + 2)^ ^(„„ 



2(2£+ 1) 



'di-l 2 



i(^-i)^<2-\ 



(B25) 



where I have defined JifJ-' 



and 



1. Angular integrals 

The calculation of the boundary integral in Sec. V B re- 
quires the evaluation of numerous integrals over a cross- 
section of the worldtube F. This subsection compiles 
these integrals. Let ~ {t, r, 6, (f>) be a point lying out- 
side the worldtube, and let a;"' = (t',.^, 6', 4>') be a point 
on the worldtube; here r and ^ denote the Fermi coordi- 
nate distances to and . Let n"' and n° be Cartesian 
unit vectors defined by the angles {0,(j)) and {9',4>'), re- 

.^2 _ 2r^nan''' is 



the leading-order flat-spacetime luminosity distance be- 
tween and x" . For brevity's sake, I introduce the 
notation (/) = J f{9',(f)')dVL' to denote an average 
over the primed angles, where dO! = sm 9' dO' d4)' . 

I group the integrals according to the power of ^ that 
appears in their integrand: 



spectively. The quantity ^ = 



1 



1 

r 

Sin'' 



(B26) 

(B27) 
(B28) 
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r2(r 



r3(r2 -^2) 
r4(r2 -^2)' 
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(B31) 
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(B34) 
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(B37) 

(B38) 
(B39) 
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APPENDIX C: EXPANSIONS OF BITENSORS 
NEAR THE WORLDTUBE 

I present here the expansions of various important 
bitensors of the form Tix^ x'), where x is a point in the ex- 
terior of the worldtube T, and x' is a nearby point on the 
worldtubc. The expansions are based on standard meth- 
ods of covariant expansion reviewed in Rcf. [29] . To keep 
track of the orders of the expansions, I use the quantity 
A ^ r ~ ^ ~ At, where r is the radial coordinate at x, 
!M is the radius of F, and At is the proper-time difference 
between x' and x. For brevity, I introduce the shorthand 
notation a;"^ = a' = a{x',x), a = a{x,x), and 

a" = (t{x,x"), where x and x" are points on the world- 
hne 7. The relationship between the various points will 
eventually be identified with that depicted in Fig. 4. 

For completeness, the expansions in this section al- 
low for an arbitrary, non- vacuum background and an ar- 
bitrarily accelerating worldlinc. In the calculations in 
Sec. VB, both the acceleration and the Ricci tensor can 
be set to zero. 



1. General expansions 

Consider a bitensor Ta'{x,x'). We can expand this in 
a sequence of steps: First, we expand along a geodesic 
that connects the point x' on the worldtube to a point 
X = 7(i') on the worldline; this is an expansion in powers 
of 1%. Next, we expand up along the worldline, about a 
point x" = 7(i); this is an expansion in powers of the 
proper-time difference At = t — t'. These two expansions 
leave us with bitensors that depend on x and x" . The 
final step in our procedure is a near-coincidence expan- 
sion of these bitensors; this last step is an expansion in 
powers of r. This procedure does not rely on any par- 
ticular relationship between x' and x or between x and 
x" . It becomes a coordinate expansion by fixing these 
relationships: for example, by connecting x to x" with a 
geodesic perpendicular to the worldline (and doing like- 
wise for x' and a;), the covariant expansion becomes an 
expansion in Fermi coordinates. 

We first hold x fixed and expand the x'-dependence 
about x: 



Ta'{x,x') =g2'Yl 



k>0 



til 

k\ 



(Cl) 

where the reader is reminded that the parallel propaga- 
tor is given by g^, = e^e^, . Next, still holding x fixed, 
we expand each of the bitensors Ta;^^...^f.{x, x) about x" . 
Since wc do not possess a conventient expression for the 
parallel propagator between x and x", we perform this 
expansion along the worldline. We do this by express- 
ing Tg in terms of its tetrad components, converting to 
tetrad components via the relationship Tg ~ Tj{t')eJg^, 
and then expanding the tetrad components in powers of 
the proper time interval At. The time-derivatives along 
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the worldline are evaluated covariantly by rcexpressing 
the tetrad components m terms of the coordinate basis, 
leading to 



3" h 

a"'" = —e)^ X , and identifying t with Fermi time. Anal- 
ogous identities would generate an expansion in terms of 
retarded coordinates {tret,x'^et) advanced coordinates 

{tadv, ^adv) 

The generalization of this procedure to tensors of other 
ranks is obvious. 



n>0 

dt" •-Ti—yki ^Ji 



(C2) 



This can be expressed in terms of covariant derivatives 
of T and combinations of tetrad and acceleration vectors 
by using the identity 



dt" ) V " -''"■■■'^fc ^ ■-'i 



(C3) 



where the derivative of the tetrad is given by ^e" = 
{u°'af3 — a"uis)ej. The indexed tensor A^^ '"^i ihj) is con- 
structed from the four- velocity and its derivatives. 
Explicitly, 



A(0,0) = 1, 



(C4) 



+ l)u^'^, (C5) 

where < j < and A{i, j) = for j < and j > i. 

Finally, the bitensors Tai'-^'^' ,,,^'^s{' ...s" [x, x") can be ex- 
panded about x" using the near-coincidence expansion 



m>0 



cr"Mi . . . cr"A'„ 



(C6) 



where tV^ ^ „ ,, (x") is defined by the recurrence rela- 
tion [49] 

t^...)(x")=£m,r(...)(a;,a;"), (C7) 
K..Wl..V:S^") = lim^Tf...),^....^. (C8) 



X — >x' 
m— 1 



V {t"\ 
^ \ f t^i ■■■t^'e''t^'e'+i---t^'m^ '' 



1=0 

Substituting Eq. (C6) into Eq. (C3), Eq. (C3) into 
Eq. (C2), and Eq. (C2) into Eq. (CI), we arrive at an 
expression for Ta'(x,a;') in terms of tensors at x" and 
the small expansion quantities ct'". At, and ct"" . This 
procedure is valid in any coordinate system. It can be 
made into a coordinate expansion in terms of Fermi co- 
ordinates (i, a;") by using the identities tr'^ = —e^x'^ and 



2. Expansions of At, (t^i{x,x'), and related 
quantities 

From this point on, I restrict myself to the case in 
which X and x' are connected by a unique null geodesic. 
In other words, x' lies on the surface S in Fig. 4. I begin 
by expanding u(x,x')] since a;' G 5, we have (t(x,x') = 0. 
We can make use of this fact to find At, which is required 
for all other expansions. 

So, following the procedure outlined above, we first 
expand about x: 



la . 1- -la -in 1- =.la-l8-M 



cr(x, x') ~ a — Oao'^ -f- -^a^^a'^-a 

_ _-n 

2A^al3^S 

We next expand a about x": for example 



^a,^,,a'^a'%'^a'^ + O(A^) . (C9) 



+ 0(A^). (CIO) 



2' ' "p" 



e^iiXa, and the stan- 



Using a" 

dard near-coincidence expansion = — 

^Rfi"-y"}^"S"0'"''" o-"^" + 0{X^), and dropping terms of or- 
der a^, we arrive at the expansion 

a = - i (1 + a,{t)x- + ^RocOd{t)x^'') (At)^ 

+ i (dc(t)x^) (At)3 + 0(a5) . (Cll) 

The same procedure yields 

^a = -Xael + {el + aa{t)x''el + ii?oa06(t)a^"^e° 
+ ii?ca06(0a;°'4)At 



hA\^-ob{t)x''ei + aa{t)e%) {Mf 



(C12) 



+ \a,{t)el{Mf + 0{\^) 

^aP = 9a0 - ^RlaJb{t)x"'''eJ^ej 

+ ^Raijb{t)x''el^ej^At 

- ii?aobo(i)e|e|(At)2 + 0(a3) (C13) 
^apj = -jRKija{i)x°'e\g^e'^^^e^ 

- lRKiJo{t)el^ejf^At + ©(A^) (C14) 

2 D „J „K„L 



'aP'yS 



iRKijL{t)e[^ejf^e'l + 0{\). 



(C15) 
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Substituting these expansions into (C9) and setting the 
result equal to zero, we get 



= cr(x, x') 

+ |i?Oahc(t)(x"^'x^ + a:'V^)At 



cd^/ab 



1 + acit)ix'= + x"=) 



(C16) 



+ ^Ro.odmx'" + x^x" + x'^")] (At) 
+ lac{tW + 2x"'){Mf + O(A^) . 

We next expand At as 

At = A (Ato + AAti + A^Ata) + 0(A*) 



(C17) 



Substituting this into (C16) and solving order by order, 
we find 

Ato = ^ 



Ati 

At2 



: -^i^ac{t){x' + x'') 

-lt^O^RacM{t)x'''x"''' (C18) 

+ ii?0a6c(i)(x''V'-2;'°=x^). 



where 



Other useful expansions are 



afj,'{x,x')n'-' 



Hi: — XaU 

- i^aa(t)n"^(^-2^i) 
-l^RoaOtrnx'^ + 2x"')n"' 



(C24) 



+ l^Roacb{t){x'''' + 2x"'x')n'^ 
-~l.^Racbdit)x''''n'''' + 0{\^) 



and 



+ ^Roabcmx" ~ x"')x'n"= 

+ 0{X^). (C25) 



3. Expansions of Green's function 

I present here the expansion of part of the Green's func- 
tion for the case in which x, a;', x, and x" lie within one 
another's convex normal neighbourhood. By following 
the same procedure and making use of standard near- 
coincidence expansions given in Ref. [29], one finds the 
following expansion for the direct part of the Green's 
function: 



^ = +^2- 2r^n''n'^ (C19) 

is the flat-spacetimc luminosity distance between x and 
x'. 

Using this result for At, we can now find an explicit 
expansion for any bitensor at x and x' . In particular, 

^ = a^i ^^j- can be expanded as 

^= A(^ + A^i + A^^^j) 
where ^ is given above, and 
i-i = i^ac(t)(a;"4-x"=). 



O(A^), 



^2 — — Racbd{i)' 
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Note that = — Ati, which is what we would expect in 
flat spacetinie. 

The time-derivative of ^ can similarly be expanded to 
find 



9t'^= ^ + A^i + A^^ + 0(A^ 



(C22) 



where 



i^o = -1 

i^ = -aamx'^ +X'') 



(C23) 
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And the expansion of its covariant derivative is given by 



'^ap -S' — ^I 6j e^, 



il 



7 /,/ 

</3 if 



^^Q/3^'^if 



+ H^/W^ + 0(A^) 



(C33) 
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where 



<(7 IJ 



i(8 IJ 



it 



9 /J 

■a/3 



_ffc 



-e(ae^)^L''/«W - ^e[„e^)i?KoW, (C35) 



e[A)RL' Kc{t) + ieLe;5)i?Kc(i) 



(C34) 
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(C36) 



APPENDIX D: EXPANSIONS APPEARING IN 
THE SECOND-ORDER WAVE EQUATION 



where the second superscript index in parentheses de- 
notes the power of r. Making use of the expansion of 
h ^ , obtained by setting the acceleration to zero in the 
results for lij^ found in Sec. IV B, one finds 



(D3) 



I present here various expansions used in solving the 
second-order Einstein equation in Sec. IV C. 

The part of the Ricci tensor that is quadratic in the 
perturbation h is written as 



and 



(Dl) 



I require an expansion of 6^R'^^j^[)i '] in powers of the 
Fermi radial coordinate where for a function /, 
(5^i?|^°^[/] consists of S'^Rap[f] with the acceleration 
set to zero. Explicitly, I require the coefficients in the 
expansion 
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I require an analogous expansion of E^^l 



where E^^p[f] is defined for any / by setting the 



acceleration to zero in Eap[f]. The coefficients of the l/r^ and l/r^ terms in this expansion can be found in Sec. IV C; 
the coefficient of l/r"^ will be given here. For compactness, I define this coefficient to be Eap- The tt-componcnt of 
this quantity is given by 



Ett = 2dlM,n' + \S'B,,tI 
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The ta-component is given by 



This can be decomposed into irreducible STF pieces via the identities 



ikl 
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e^aiM(,B,k) + \5aiMPB',,e 



which follow from Eqs. (B3) and (B7), and which lead to 
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+ 4eai' {S^jSki) + '^M(^Bki) ) 
The a5-component is given by 

Eab = fm^£^Jnab'' + fm^Sab - Sab [{^d^ + IS^B.j + f AP£„) + (fM^.S^k) - f ^(^S,,)) n'^" + mm^£,,n'^] 
+ ^Mi^,£,^,n' + ^M'£,(,nfe) + fm^S^^^fi,^ ' + 16A.f,£,(,ri,)''^' - f + ± {10S,Bab + 27M,£ab) 
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Again, this can be decomposed, using the identities 
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which lead to 



Eab = -2Sab + IS^B,, + lAPS^j) + fm^£,,n'^ + {fM^.S^k) - P^.B^k)) n^'"] 
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